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PREFACE 


The purpose of this contract was to perform an investigation of a 
hearingless helicopter rotor concept having a composite primary structure. 

The results of this investigation are presented in two separate reports. 

The one herein contains the results of graphite/epoxy fatigue tests, wind 
tunnel experiments, correlation studies, and a preliminary design of a 
full scale helicopter rotor. The companion report, NASA CR-2638 presents 
the mathematical derivation of and program user's manual for the aeroelastic 
analysis used to generate most of the analytical results presented herein. 

The project was managed hy Mr. M. C. Cheney. The composite material 
investigation was performed hy Mr. R. C. Novak. The detailed design and 
fabrication of the model rotor was performed hy Mr. Cheney and Mr. E. D. Bell- 
inger, and the experimental data were acquired hy Messrs. Bellinger and 
B. W. Goepner. Development of the aeroelastic analysis, generation of the 
correlation results, and the aeroelastic analysis calculations of the full- 
scale design were performed hy Dr. R. L. Bielawa. Design of the full-scale 
hlade configurations was performed hy Mr. J. A. Longohardi of the Sikorsky 
Aircraft Division with assistance from Dr. L. E. Greenwald of UTRC. Assis- 
tance in the management of the experimental test program was provided hy 
Mr. A. J. Landgrehe. 

Note that all original data was recorded in English units and suh^iequently 
converted to metric units. A table is presented to facilitate conversion of 
data hack to English units should it he desired. 
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INVESTIGATION OF A BEARINGLESS HELICOPTER ROTOR 
CONCEPT HAVING A COMPOSITE PRIMARY STRUCTURE 

Richard L. Bielawa, 

Marvin C. Cheney, Jr., and 
Richard C. Novak 

United Technologies Research Center 
SUMMARY 


Experimental and analytical investigations were conducted at the United 
Technologies Research Center (UTRC) to evaluate a hearingless helicopter rotor 
concept (CBR) made possible through the use of the speeieQ.ized nonisotropic 
properties of composite materials. The investigation was focused on four 
principal areas which were expected to answer important questions regarding 
the feasibility of this concept. First, an examination of material properties 
was made to establish moduli, ultimate strength, and fatigue characteristics 
of unidirectional graphite/epoxy, the composite material selected for this 
application. The results confirmed the high bending modulus and strengths and 
low shear modulus expected of this material, and demonstrated fatigue proper- 
ties in torsion which make this material ideally suited for the CBR application. 
Second, a dynamically scaled model was fabricated and tested in the UTRC low 
speed wind tunnel to explore the aeroelastic characteristics of the CBR and 
to explore various concepts relative to the method of blade pitch control. 

Two basic control configurations were tested, one in which pitch-flap coupling 
co\ild occur (cantilever torque tube design) and another which eliminated 
all coupling (pinned-pinned torque tube design). It was found that both 
systems could be operated sucessfully at simulated speeds of 180 knots; 
however, the configuration with coupling present revealed a potential for 
undesirable aeroelstic responses. The uncoupled configuration behaved 
generally as a conventional hingeless rotor and was stable for all conditions 
tested. Another decoupled rotor configuration (snubber torque tube design) 
was extensively analyzed, but was not tested due to scaling difficulties. 

Third, a preliminary design of a full-scale rotor was conducted to establish 
full-scale requirements and to show that no significant problems would be 
encountered in the manufacture of a CBR blade. Finally, an aeroelastic 
computer program was developed to sim^llate the unique CBR blade structtire and 
operating environment. This analysis was then employed to examine the sta- 
bility characteristics of the full-scale snubber torque tube design under 
various flight conditions and under conditions where material properties were 
degraded due to possible flaws in fabrication or partial failures in operation. 
This study revealed no dynamic problems for the conditions investigated, has 
confirmed that significant reductions in rotor system complexity and weight 
can be achieved with the CBR, and further, that graphite/epoxy is ideally 
suited for the special requirements of the CBR. 



A description of the aeroelastic analysis developed under this contract 
along with detailed procedures for its use is presented separately in NASA 
CR-2638 (Eef. 1 ). 
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INTRODUCTION 


Fully articulated rotors have successfully powered rotary wing aircraft 
since their inception and have brought the helicopter to its present unique 
position in the field of air transportation. However, high cost and mainte- 
nance requirements, primarily related to complicated hub mechanisms, have 
restricted the growth of the helicopter. Hingeless rotors sought to reduce 
these requirements by simplifying the hub through the elimination of the 
flapping and lagging hinges; however, current designs (Refs. 2 through 4) 
still have not significantly reduced weight and cost since feathering bearings 
have been used to provide pitch control and heavy structure has been necessary 
to react centrifugal loads. 

The study of a new rotor concept designed to eliminate all hinges and 
bearings through a specialized use of fiber reinforced composite materials, 
was initiated at the United Technologies Research Center (UTRC) in I968. This 
concept utilized a finite length of unidirectional composite material, with 
its attendant low shear modulus, in place of the feathering bearings of con- 
ventional hingeless rotors. This "flexbeam" member is carried through the 
center region to retain an opposing blade, thus eliminating individual blade 
retention hardware and resulting in a simple, low weight, low profile hub 
requiring little or no maintenance. The encouraging results obtained from the 
early work at UTRC led to a research contract jointly sponsored by the NASA^" 
Langley Research Center and the US Army Air Mobility Research and Development 
Laboratory. It is the results of this contract that are reported herein. 

The objective of this research was to assess the general feasibility of this 
rotor system, called the Composite Bearingless Rotor (CBR) through a series 
of experimental and analytical investigations. It was the p\irpose of these 
investigations to explore all aspects of the rotor system, from the viewpoints 
of basic material properties, stability characteristics, wind timnel verifica- 
tion of rotor performance and dynamics, and full-scale design requirements. 
There were fundamental questions regarding the durability of unidirectional 
composite materials when utilized in the unconventional manner required of 
this design, that is, with imposed oscillatory motions about the weakened shear 
axis where the fibers are not loaded in their respective axial directions. 
Another potential problem area which required study was the effect a time 
varying twisted beam mi^t have on the blade response. In effect, the inboard 
blade structure varies with collective and cyclic pitch potentially influencing 
the in-plane and out-of -plane stiffness characteristics. The problem of in- 
putting control deflections by means of a torque tube at a finite distance out- 
board of the hub also had to be considered, as well as pitch-flap and pitch-lag 
coupling phenomena which may result from motions of this control torque tube. 
Finally, the problem of applying the technology developed and the experience 



gained in the model scale investigations to full-scale rotors had to he 
addressed. The task of eventually reducing these results to practice is, in 
the final analysis, the ultimate test of the success of a new concept. 

A review of the CBE concept and a svmnaary of a portion of the resiolts 
obtained during the course of the contract were presented at the American 
Helicopter Society National Forum in 1972 (Ref. 5). 
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LIST OF SYMBOLS 


aerodynamic section lift curve slope, (deg)“^ 

longitudinal and. lateral cyclic pitch angles, respectively, 
applied at the swash-plate, deg 

nimiber of blades 

blade or flexbeam section chord, m 

airfoil section lift, drag and pitching moment about 
quarter chord aerodynamic coefficients, respectively 

airfoil section drag aerodynamic coefficient at angle-of- 
attack for zero lift 

rotor drag coefficient, D/pTrR^(f3R)^ 

rotor lift coefficient, L/pTTR^(fiR)^ 

rotor propulsive force coefficient, PF/pnR^(nR)^ 

pitching moment coefficient, My/pTTR^(nR)^ 

rotor torque coefficient, 0/pnR^(nR)^ 

rolling moment coefficient, Mjj/pnR^(OR)^ 

rotor thrust coefficient, T/pnR^(fJR)^ 

lateral force coefficient, y/pttR^(CB)^ 

•3 p 

rotor drag, N; or plate bending stiffness, Et /12(l-u ), 

Nm 

radial coordinate of the blade offset point 

p 

Young’s (bending) modulus of elasticity, N/m^^ 

blade bending stiffnesses for edgewise and flatwise 
flexing, respectively, 



LIST OF SBIBOLS (Cont'd) 


aei 


""PR 

G 

GJ 

l/c 


L 

Me> Mp 


Mx» My 

^0 

PF 

nP 

Q 

r 

R 

R 

[S] 


difference in blade bending stiffnesses, EI^-Elj, ,Nm^ 

torsion shear stress, MN/m^ 

push rod load, positive in compression, N 

shear modulus of elasticity, N/m^ 

2 

torsional stiffness, Nm 
section modulus, m3 

2 

torsional inertia distribution, kg m /m 

radius of gyration of tension load carrying portion of 
section, m 

rotor lift, N 

blade edgewise and flatwise bending moments, respectively, 
positive forward and upper fibers in compression, Km 

rotor aerodynamic roll and pitch moments, positive port and 
nose edges up, respectively, Nm 

torsion moment, positive for radially increasing resultant 
twist rate. Km 

rotor propulsive force, positive in direction of flight, N 
n times rotor frequency, n 
rotor torque, Nm 

blade (nondimens ional) radial station, x^/r 
rotor radius, m 

ratio of full scale rotor radius to model rotor radius 
aerodynamic static stability matrix 
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LIST OF SYMB0I6 (Continued) 
blade or flexbeam section thickness, m 
rotor thrust, K 
rotor aerodynamic trim matrix 
local aerodynamic section air velocity, m/sec 
forward fli^t speed of rotor, m/sec 
blade mass distribution, kg/m 
spanwise location, m 

blade (nondimensional) inplane and out-of-plane elastic 
deformations , respectively 

rotor lateral force, positive to starboard, K 

airfoil section angle of attack, deg 

rotor shaft angle, positive is nose-up deflection of 
plane of rotation relative to tiinnel air velocity, deg 

fiber alignment angle for composite materials, deg 
(see figure 1) 

uncoupled vibratory mode shapes for i^^ flatwise bending, 
kth edgewise bending and torsion modes, respectively 

effective flappixig hinge orientation angle producing 
pitch-flap coupling, deg 

blade equivalent mass ratio, pacR/j^ wCYw^) dr 

ratio of damping level to critical value 

blade collective control angle experimentally measured 
at swash-plate or other remote location, deg 

automatic blade pitch change per xmlt deflection in the 
ith flatwise bending mode 
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LIST OF SYMBOLS (Continued) 


blade collective angle as measured at 0.75 station, deg 

blade linear twist angle defined by twist rate of 
aerodynamic portion, positive for root pitch less than 
tip pitch, deg 


^®.75R’ 
^s 
X 9 Xq 


perturbations in collective, longitudinal and lateral 
cyclic, and rotor shaft angles, respectively, deg 

uniform rotor inflow ratio, also eigenvalue, as appropriate 


Xt 9 X-i 

■*-C J-s 


first harmonic cosine and sine components of radially 
linear rotor inflow, as defined by Glauert 




advance ratio, Vcosofg/OR 


V 

P 

a 




T 

(I) 




f 




Poisson's ratio 

O 

air density, kg/m'^ 

rotor solidity, bc/rrRj real part of eigenvalue, or modal 
equivalent characteristic exponent 

blade edgewise and flatwise bending stresses, respectively, 
positive for leading edge and upper surface in compression, 
m/rn^ 

torsional shear stress, N/m^ 

blade azimuth angle measured from downstream position, deg 

blade nattiral frequency (for uncoupled or coupled mode); 
imaginary part of eigenvalue; or identified modal frequency, 
as appropriate, Hz 

blade (nondimens ional) natural frequencies for ith flatwise 
bending, kth edgewise bending and jth torsion uncoupled 
nattiral modes, respectively 
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r 


Cl 

n 


Subscripts 
( )e 

( )fb 
( )ps 

( )f 

( )i 
^ ^nom 


^ ^subcr. 

^ ^supercr. 

( )trim 
( )tt 

Superscripts 

(“") 

(*) 

( )’ 


LIST OF SYMBOLS (Continued) 
rotor stational speed, rad/sec 

ratio of fvill scale rotor rotational speed to model rotor 
rotational speed 

due to elastic effects 
for flexbeam 
full scale 
final value 
initial value 

model scale 

for nominal configuration 

for subcritical design configuration (i.e. > 1) 

for supercritical design configuration (i.e. s 1) 
trimmed conditions 
for torque tube 

nondimens ionlizat ion by combinations of R and Cl 
differentiation with respect to 
differentiation with respect to r 
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CONCEPT DESCRIPTION 


The Composite Bearingless Rotor (CBR) evolved as a resiilt of the unique 
application of fiber reinforced composite materials. The nonisotropic nature 
of unidirectional composites, generally considered an undesirable feature, 
yields a structure which was decidedly softer in torsion than in a direction 
which loads the fibers axially. By en5>loying a finite length of this material, 
torsional deflections can be achieved to satisfy the pitch angle change require- 
ments of a helicopter rotor while providing sufficient stiffness in the in- 
plane and out-of -plane directions to satisfy hingeless rotor frequency require- 
ments. 

A comparison of the modulus characteristics of fiber-reinforced composites 
and metals is given in Pig. 1 where Yoving's modulus, E, is plotted versus shear 
modulus, G. The variations, or trade-offs, of E and G for composite materials 
can be achieved by adjustments in the alignment of the fibers relative to the 
spanwise axis of the structure denoted by P in Fig. 1. The shaded area rep- 
resents a boron composite where the left end corresponds to material with 
unidirectional fibers (P = O) producing high bending modulus and low shear 
modulus, and the right end of the shaded area corresponds to a material with 
all fibers aligned ±45 deg giving maximm shear modulus. The E and G values 
between these extremeties are achieved by alignment angles between 0 deg and 
45 deg or by combinations of unidirectional and cross plies. Values are 
shown for several metals which demonstrate the increasing trend of shear 
modulus with Young's modulus. This is typical of isotropic materials which 
have a fairly constant value of E/G of approximately 2.5 as shown by the 
straight line in Fig. 1. Composite materials can have the opposite trend - 
decreasing G with increasing E which is the key to the use of these materials 
in the bearingless rotor concept. 

A schematic of the basic components of the CBR are shown in Fig. 2. The 
blade spars consist of unidirectional composite material as shown in Fig. 2a, 
where, for a 4-bladed rotor, two single lengths are used, one for each set of 
opposing blades. The 4-bladed system simply consists of intersecting the two 
spars as shown in the figure. The dimensions of the spar are dictated by the 
edgewise and flatwise frequency requirements and the required torsional prop- 
erties over the inboard, or flexbeam, portion of the blade. The aerodynamic 
cover is assembled from approximately the 25 percent radial position to the 
tip. This cover would normally be constructed of composite material with 45 
deg fiber alignment to give maximum torsional stiffness. The cover is then 
bonded to the spar through a foam filler or honeycomb resulting in a torsion- 
ally stiff blade over the outer 75 percent as shown in Fig. 2b. The inner 
25 percent, the flexbeam, remains soft in torsion, and it is this portion of 
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the blade that is elastically twisted to achieve collective and cyclic pitch 
control. The remainder of the rotor assembly consists of the control torque 
tube which is cantilivered inboard from the 25 percent station. Fig. 2c, and 
the two plates which make up the h\ib. Fig. 2d. Because of the carry-thro\igh 
feature of the blades, minimal h^ib retention hardware is required resulting 
in an extremely simple hub design. 

There can be several variations to the CBR concept which are not depicted 
in the sketches in Fig. 2. For example, it may not be necessary to extend 
the unidirectional spar over the full blade span. A conventional blade design 
could be employed and joined to the outboard end of the flexbeam. Also, 
dynamic considerations may require more or less fixity on either the inboard 
or outboard end of the torque tube than is inherent in the cantilevered design 
shown for concept illustration purposes in Fig. 2. For example, the inboard 
end of such a cantilevered torque tube could be additionally attached to the 
flexbeam by means of a "snubber" which woTild allow relative feathering motion 
but suppress relative translational motion. Various such alternative torque 
tube configxirations are described in detail in subsequent portions of this 
study. 

Controls to the rotor are applied through a swash plate in a conventional 
manner. The dynamic response characteristics are generally similar to hinge- 
less rotors operating today, particvilarly for the snubber design where nearly 
all blade pitch coupling is eliminated. These response characteristics are 
almost entirely dependent upon the fundamental bending frequencies which can 
be tailored to any desired values without affecting the basic CBR concept. 
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COMPOSITE MATERIAL INVESTIGATION 


/ 


As discussed previously, the bearingless composite rotor concept is 
heavily dependent on the unique elastic properties of the materials as well as 
their fatigue characteristics. As a result, high priority was placed on a com- 
posite material investigation to define modulus characteristics, measiire static 
strengths, and determine fatigue properties. 

The fatigue performance of advanced resin matrix composites is generally 
considered to be a plus when the materials are compared with metals. The bulk 
of the supporting evidence for this conclusion is based on testing performed 
vtnder conditions in which the reinforcing fibers carried the applied load 
(Refs. 6 through 8). Thus, since the fibers behave in an elastic manner, it 
is not suprising that the fatigue performance of the composites has been foimd 
to be good. However, the bearingless rotor concept relies on the fatigue per- 
formance of the composite structure under combined bending and torsion loading 
conditions. Under these circumstances the composite response is not necessarily 
elastic over all load ranges, and the matrix resin strength or the resin-fiber 
interfacial bond strength becomes very important to the performance of the 
system. 

The following sections describe the materials and test procedures used in 
an experimental program to study the pertinent static and fatigue properties 
of the composites, and a brief discussion of the test results. 


Experimental Approach 

All specimens consisted of HTS (high tensile strength) graphite fiber 
in a matrix of Fiberite X-05 epoxy. Preimpregnated tape (prepreg) was 
purchased at various intervals in the program to provide information 
on the repeatability of material properties from lot to lot. The speci- 
mens were fabricated from 12 plies of prepreg and cured at 325^F under 
.55 MN/m^ (80 psi). The finished specimen sheets were 2.5 nim (O.IO in.) 
thick with a fiber content of approximately 60 percent. The fatigue test 
specimens were cut to either 2.5^ cm or 1.27 cm (one inch or one-half inch) 
widths and approximately fifteen cm (six inches) in length* The tests consisted 
of static flexure and interlaminar shear, flexure and torsion fatigue, and 
combined flexure/torsion fatigue which simulated combined blade flatwise, 
edgewise, and torsional deflections. Static U-point flex\ire tests were 
conducted at a span-to-depth ratio of 32:1, Interlaminar shear strength 
was measured in a 3-point bend test at a span-to-depth ratio of ^:1. The 
fatigue tests were conducted using a VISHAY variable throw cam driven by a 
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1120 watt (1.5 hp) electric motor at 3600 rpm as shown in Fig. 3. The specimens 
were run to lo"^ cycles, and the post-test moduli were compared to the pre- 
test moduli to determine possible material degradation. 


Results and Discussion 


Static Properties - The results of the static property measurements are 
presented in Table II. With the exception of the flexural strength of the Lot 
lD-52 IB, the data exhibited fairly low scatter and fell within the e^qpected 
range. The apparent low strengths of that lot were due to a low filament 
volume fraction (50 percent) in the composites which were tested. This was also 
reflected in the somewhat lower bending mod\ilus for the material. 

The pertinent points bro\ight out bF the static data are the verification 
of the required high E/G ratio (~25) and the demonstration of static shear 
and bending strengths which are significantly higher than the design stresses 
for a typical bearingless rotor application. 

Fatigue Properties - The res\ilts of the pure bending fatigue and pure 
torsion fatigue are graphically presented in Figs. 4 and 5» respectively. 

As mentioned previously the performance of the materials was measured by 
reductions in the shear stiffness, G, and the bending stiffness, E, 
because the specimens did not suffer complete fracture as metals do under 
fatigue loading. Rather, the failure mode was in the form of cracks 
throiagh the resin matrix which gradually reduced the section modulus of the 
specimen. An example of the cracking which was observed is shown in Fig. 6 
which is a photomicrograph of a specimen subjected to a vibratory torsional 
stress of ± 52.7 MN/m^ (?650 psi) for 5 x 10° cycles. The graphite 
filaments appear as light circles which are surrounded by the epoxy matrix. 

The small crack can be seen to propagate from interface to interface 
through the matrix. There is no evidence of broken filaments. 

Since no precipitous stiffness reduction was evident from these results 
it was difficult to establish a failure criterion. Accordingly, the stress 
corresponding to a nominal 10 percent reduction in modulus was tentatively 
selected as the endurance limit. Considerably more fatigue testing wotild 
be necessary before a firm endturance limit covild be established. The vibratory 
stresses corresponding to the 10 percent criterion are estimated from Figs. 4 
and 5 at 758 MN/m^ (110,000 psi) for bending and U8 MN/m^ (7000 psi) for 
shear. 
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Because the vibratory stresses are present in combined form in the 
bearingless rotor, a series of combined bending/torsion fatigue tests was 
conducted to determine if the material responded in a predictable manner 
based on the pure bending and torsion results. In addition, some tests were 
performed with a steady axial load to simulate the centrifugal force present 
in an actual blade. 

The combined stress state was produced by using a setup similar to that 
shown in Fig. 3 except the specimen was positioned off the center of rota- 
tion of the fixture, producing a combination of bending and twisting. Steady 
tension was initially provided by dead weight loading of a cable attached to 
the end of the specimen by friction clamps. Later this was replaced by a 
hydraulic cylinder. Several testing trials were conducted before satisfactory 
results were achieved. 

The results of these tests are summarized in Table III. In general, the 
data indicate that the effects of combined loading were slightly more severe 
than would be predicted based on the simple loading results discussed pre- 
viously. This may have been related to test effects such as indicated for 
specimen KCl- 6 . Nevertheless, with the exception of specimen KCU-2, the 
specimens were capable of sustaining a total bending stress of approximately 
669 MN/m^ ( 9 T 5 OOO psi) for lo"^ cycles with 10 percent reduction in modulus. 
This compares favorably with the 758 MN/m^ (110,000 psi.) endurance limit found 
for pure bending, given the stochastic nature of fatigue. 
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MODEL ROTOR TESTS 


Test Facility and Bq,uipment 

Wind Tunnel - The tests were conducted in the United Technologies 
Corporation 1.2?.- x 1.83-ni subsonic wind tunnel which has a cilpsed return, 
closed throat circuit and operates with a maximum test section velocity and 
Reynolds number of 44.7 m/sec and 2.89 x 10^ per m, respectively. Static 
pressure in the irregular octagonal test section is atmospheric. 


Small-Scale Rotor Test Rig - The rotor test rig employed in this program is 
shown in Fig. 7. It was designed primarily for exploratory testing of small- 
scale rotors in the 1.22- x 1.83-m subsonic wind tunnel. Power to the rotor 
is supplied by direct drive from a 40 hp, variable-speed electric motor. The 
motor and rotor force measuring system were mounted outside the tunnel above 
the test section. Rotor thrust, drag and torque were measured by means of 
strain-gaged bending beams. A 25-unit slip-ring assembly mounted on the rotor 
drive shaft provided a signal transfer point for measuring blade stresses and 
bending loads. 


Model Rotor Description 

Both 2- and 4- bladed dynamically scaled model rotors were tested in 
hover and forward flight to demonstrate the viability of the composite bearing- 
less rotor concept. The blade design parameters for the two torque tube con- 
figurations are presented in Table IV. A photograph of a 4-blade rotor in- 
stalled in the wind tunnel for forward flight test conditions is shown in 
Fig. 7 . The elements comprising the flexbeams and the pinned and cantilever 
torque tube configurations are shown in Figs. 8 and 9s respectively. 

The model rotor blades were fabricated as integral sets of opposing blade 
pairs. Each such blade pair set utilized a common primary structural bending 
element comprised of unidirectional graphite filaments, continuous from tip 
to (opposite) tip, cast in an epoxy matrix. Furthermore, each rotor blade in 
turn consisted of two main spanwise sections starting from the blade root 
offset (rotor hub). The first spanwise section comprised the fl exbeam -torque 
tube region and the second made up the outer blade region. The details of the 
flexbeam and two torque tube configurations are shown in Figs. 8 and 9- 
flexbeam, which was the inboard portion of the centrifugal and bending load 
carrying member, had zero pretwist, spanned the radial region from O.05 to 
0.33R and had a cross sectional chord to thickness ratio of approximately 7.2. 
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Outboard of the flexbeam, additional model spar components consisted of the 
graphite /epoxy material encased in aluminum 0.25^ mm thick to produce the model 
scale torsional stiffness. Lead tape 0.127 mm thick was used to scale mass 
and inertias, and balsa wood was used to form the NACA 0012 airfoil. A pro- 
tective wrap of 0.064 mm mylar was shrunk fit to the balsa to provide an 
aerodynamically clean and uniform surface. Figure 10 illustrates the 
Juxtaposition of the various elements of the outboard blade cross section. 

The blade chord was constant at 0.039 ® from the 0.33 R station to the 
tip. The blade was preconed 2 deg to alleviate the steady bending stress due 
to rotor thrust at a C^/ct = 0.08. The model blade spanwise distributions of 
wei^t, flatwise stiffness, edgewise stiffness, torsional stiffness, torsional 
Inertia and the spar radius of gyration are plotted in Figs. 11 and 12 for 
both pinned and cantilevered torque t\ibe configurations. Using the blade 
section properties shown in Figs. 11 and 12, the uncoupled model blade (vacuum) 
natural frequencies and mode shapes were calculated using the United Tech- 
nologies Corporation Uncoupled Mode Blade Watriral Frequency computer program. 
These are shown in Figs. 13 through l6 and are given in Table IV for the 
nominal tip speed conditions of 99*1 m/sec. The flatwise and edgewise bending 
modes (Figs. 15 and l6) show the typical character of a cantilevered beam 
with zero deflections and slope at the blade root offset. 

As stated above, two torque tube configurations were tested.on the CBR 
blades. The pinned-pinned torque tube configuration, shown in Fig, 8, con- 
sisted of a pitch 6irm, feathering bearing, universal Joint, 0.635 cm OD- 
0.089 cm wall thickness 303 stainless steel tubing, a flexible coupling and 
blade attachment. This torque tube allowed no pitch-flap coupling because 
the universal joint and flexible coupling transmitted only pitch moments 
to the blade and eliminated any blade bending feedback (pitch-flap coupling) 
to the pushrod. The root end of the pinned-pinned torque tube was attached 
to the hub so that the blade centrifugal load was not significantly affected 
by the torque tube. 

The cantilevered torque tube configuration shown in Fig. 9» consisted of 
a 0.953 cm OD - 0.124 cm wall thickness 6065-T6 aluminum tubing bonded to 
two 2.54 cm square 6065-T6 aliiminum plates. These plates, together with 
their connecting clamping screws, were used to effect attachment to the blade 
at the 0.31 R spanwise station. As shown in Fig. 9) 'the cantilevered torque 
tube is mounted in its trailing edge position. It could also be mounted in 
a leading edge position. Since the cantilevered torque tube to blade attach- 
ment point was fixed at the 0.3I R station, pitch- flap coupling was changeable 
only by attaching the pitch link at various of the radial points on the torque 
tube as shown in Fig. 9* 
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Scaling Factors 


To apply the model scale data presented herein to a full-scale rotor 
of the same design, an appropriate scale factor must he applied to the 
parameter being scaled. A summary of these scale factors is given in Table V. 


Test Data 

Test Conditions - Wind timnel test conditions were selected to encompass a 
range of advance ratios, and thrust coefficients typical of contemporary heli- 
copter operation. Advance ratios of 0, 0.1, 0.25, 0.35 ^-nd 0.47 were tested 
at the nominal tip speed over a range of shaft angles and collective pitches. 
Except where specifically required for derivative data, the cyclic Aj^^ and 
control angles were adjusted at each test point to approximately mall out the 
first harmonics of inboard flatwise bending stress. Hovering dataware taken 
within the tunnel test section at 107.7 percent, 100 percent, 92.3 percent 
and 76.9 percent nominal tip speed of 99-1 m/sec for blades configured with 
the pinned-pinned torque tube and with the cantilevered torque tube incorpo- 
rating various available pitch-flap couplings. In particoolar, those canti- 
levered torque tube configurations selected included push- rod attachment points 
at 0.060 R and 0.080 R for trailing edge attachment and at 0.070 R and 0 . 09 I R 
for leading edge attachment. These four push-rod attachment configiorations 
define equivalent (geometric) 6_ angles, relative to the first flatwise bending 
mode, of 48.1°, 33.2°, -4l.5° and -23.0°, respectively. A soommary of the 
forward flight test conditions is given in Table VI. 

Test Data Acquisition and Reduction - The rotor data obtained during the 
-test consisted of rotor performance (thrust, drag, and torque) and blade 
stresses (flatwise, edgewise, and torsional) . The thrust, drag and torque data, 
obtained from self-balancing potentiometers, were corrected for tares and re- 
duced to coefficient fonn {CjJa, C-q/o and Cq/ct) using the measured shaft 
angles (hub plane angles of attack) with appropriate corrections. Flexbeam 
flatwise and edgewise dynamic stress, blade outboard flatwise bending moment, 
and flexbeam torsion moment (or alternately, twist deflection) data were 
obtained from various strain gages and numerically calculated from visl- 
corder output. These strain gages were located as follows; flatwise and edge- 
wise bending gages were located near the blade root at 6 percent span, a 
flatwise bending gage was located outboard of the torque tube attachment point 
at 42 percent span and a torsion gage was located on the flexbeam at 17 per- 
cent span. The bending stress and moment calibrations were conducted by 
applying bending moments and dividing by the appropriate l/c's. Additionally, 
the directly calculable stresses and moments due to gravity were siobtracted 
out. The flexbeam torsion gage output was similarly calibrated by applying 
a torque with the torque tube in place but with the pitch link disconnected. 
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Nominally, collective dnd cyclic pitch control angles were obtained 
from the outputs of potentiometers geared to the swash-plate drive motors. 

A cross-check on the pitch control angles was afforded by taking the mean 
and first harmonic cosine and sine components of the flexb.eam torsion gage 
output. For the pinned-pinned torque tube configuration the potentiometer 
measurements of control angle correlated reasonably well with the torsion 
gage measurements; whatever discrepancies present were most probably due 
to flexibility in the flexible coupling and/or slop between the feathering 
cuff and its retention bolt (see Fig. 8). For the cantilevered torque tube 
configuration, however, the potentiometer outputs were inherently incapable 
of giving accvirate measurements of the effective control angles due to the 
substantial washout caused by flexbeam bending flexibility. Hence, the 
potentiometer measured values of control angle used to identify the various 
performance and stress data must be considered to be only "indicated" or 
"swash-plate" values and not effective control angles. 

Rotor angle of attack was obtained by first measuring the geometric 
tilt angle of the shaft relative to the tunnel test section and then sub- 
tracting that tilt angle required, at each advance ratio, to null the thrust 
when all control angles are set to zero. For advance ratios of 0.25, 0.35> 
and 0.47, the tilt angles so required were 0.3 deg, 0.7 deg and 1.0 deg, 
respectively. The source of this zero thrust angle of tilt is the blockage 
caused by the rotor rig. Thus, this angle of attack correction is only 
approximate since a finite rotor thrust would necessarily change the rotor 
drag (or alternately, the propulsive force) which woiild in turn vary the 
blockage. However, an assessment of the variation of this blockage with 
rotor thrust was beyond the capabilities of tunnel instrumentation and the 
zero thrust tilt angle corrections were used throughout. 

Accuracy of Data - The estimated accuracies with which the parameters 
determining a given test condition could be set or determined are as follows 

parameter Accuracy 


Swash plate (indicated 


Control Angle) 

± 0.2 deg 

Tip Speed 

± 0.3 m/sec 

Shaft Angle 

± 0.05 deg 

Tunnel Speed 

± 0.3 m/sec 

Air Temperature 

± 2°F 
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The static data repeatability for thrust, drag, and torque was derived 
from repeated calibrations of the strain gages made while determining the 
calibration derivatives for conversion of strain gage units (SGlTs) into 
forces or moments. The dynamic data repeatability was established by con- 
sidering the range of Cj/a, Cjj/ct and Cq/c meas\xrements observed from con- 
secutive test points. The mean range was established and related to the 
standard deviation using the procedures outlined in Ref. 9 . The res\alts 
for the accuracy of the thrust drag and torque measurements are given below; 


Thrust, Drag and Torque Accuracy 


Static Repeatability 
Thrust: ± O.O68 N 
Drag: ± O.O9I N 

Torque; ± O.OI38 Nm 


Dynamic Repeatability 

± 0.002 Cj^CT 

± 0.004 Cp/cj 
± 0.003 Cq/<j 


The estimated accuracies of the dynamic data which were acquired on 
the recording oscillograph are the following; 


parameter 

Flatwise Bending Stress (r.= O.06 R) 
Edgewise Bending Stress (r = O.06 R) 
Flatwise Bending Moment (r ■ 0.42 R) 
Flexbeam Torsion (r = 0.I7 R) 


Accxiracy 

± 0.l4 m/xS 

± o.i4 m/m 

± 0.00043 Nm 

( ± 0.00043 Nm 
(± 0.1 deg 


Experimental Results 


The experimental results presented herein consist of hover performance 
(Figs. 17 through 22), forward flight performance (Figs. 23 through 33 )> 
vibratory stresses (Figs. 34 through 42), and control derivatives (Table VII 
and Figs. 43 through 47). Within each category the res\alts for the pinned- 
pinned torque tube are presented first, followed by those for the cantilevered 
torque tube with various of the four tested push-rod attachment locations. 
Since the pinned-pinned torque tube configuration is inherently free of any 
geometrically induced pitch- flap coupling effects, the results obtained for 
this configuration provide a good standard for evaluating the restilts obtained 
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with the cantilevered torque tube conf igxirations . All performance results 
are presented in the form of nondimens ional force and torque coefficients 
divided by solidity. The vibratory stress (and bending moment) results are 
presented in the form of 1/2 peak-to-peak (PTP) values. The "performance" 
control derivative results are presented in the form of tables of nmerically 
c8^.culated derivatives. The ";stress" control derivative results are presented 
in the form of (azimuthal) time-histories of blade stress and their perturba- 
tions due to perturbations in the control angles. 

Hover Performance - Figures 17 through 21 show for each torque tube con- 
figuration the parametric variations in CT/o’^nd Cq/t with collective pitch 
angle for each of four (4) tip speeds. Figure 17j which presents the resiilts 
for the pinned-pinned torque tube, shows typical rotor characteristics for 
variations in both collective angle and tip speed. Since the pinned-pinned 
torque tube is inherently devoid of any pitch-flap coupling, these resiilts 
provide a convenient basis for interpreting the results for the various canti- 
levered torque tube configurations which exhibit substantial pitch-flap 
coupling. Some scatter is seen in the torque coefficients at zero pitch and 
at the higher pitch settings. Variations in the effective pitch-flap coupling 
which the various push-rod attachment locations on the cantilevered torque 
tube configurations produce are amply shown in Figs. l8 through 21. The effec- 
tive pitch-flap coupling for a hingeless rotor must necessarily be determined 
from the geometry of the elastically deflected blade. The variations in per- 
formance characteristics with tip speed for the cantilevered torque tube con- 
figurations demonstrate that this geometry is significantly influenced by the 
blade's centrifugal force field. 

For those cantilevered torque tube configurations with trailing edge 
attachment (Figs. l8 and 20), the effects of control washout, due to flexbeam 
bending flexibility and negative (flap up-pitch down) pitch-flap coupling, 
are seen in the consistently higher collective angles required to obtain any 
given performance point. However, for those configurations with leading-edge 
attachment (Figs. 19 and 21) the effects of positive pitch-flap coupling when 
combined with the washout characteristics of flexbeam bending flexibility, are 
seen to produce initially a control washout at low collective angles, which 
diminishes at the higher angles. This can be seen by comparing thrust levels 
at low collective pitch settings on Figs. 19 and 21 with those on Fig. 17 for 
the pinned-pinned case. The latter case is seen to produce higher C,j,/ct for 
the same 0 at low values of 0, but at the higher 0's the C^/a values are 
similar. This trend is also seen in Fig. 22 where the pitch-flap coupling 
characteristics of the cantilevered torque tube on hovering performance are 
summarized. For either type of cantilivered torque tube attachment (leading- 
or trailing-edge ) there is a washout due to flexbeam bending flexibility. 

That is, any finite push-rod load will not only twist the flexbeam, but bend 
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it as well, due to the relatively long load path for the push-rod shear load. 
ThTis, for unit push-rod load, less of the total strain energy must necessarily 
be stored in torsion with less resulting flexbeam torsional deflection. Added 
to this initial washout due to flexbeam bending is the automatic pitch change 
due to the geometry of the elastically deflected blade (pitch-flap coTjpling), 
idiich is either positive or negative as the push-rod attachment point is ahead 
or aft of the feathering axis, respectively. For those configurations tested, 
the cantilevered torque tribe generally produced a net washout of the inputted 
swash-plate collective angle relative to the behavior of the pinned-pinned 
torque tribe. 

Forward Fli^t Performance - a simixar manner to the hovering per- 
formance resrilts, the forward flight performance results are shown in Figs. 

23 through 33 • For • each torque tube configuration Cj/(r, Oj^fcrt and C^/cr are 
shown for parametric .variations in rotor angle of attack, collective angles, 
and advance ratio. Figures 23, 24, and 25 present the performance results for 
the pinned-pinned torque tube at advance ratios of .25, .35, and .47, respec- 
tively. Dicluded in these figures are performance results for another model 
rotor similar in size, geometry, and twist and tested in the same facility 
at similar test conditions. This model rotor, however, was of semi-articulated 
type (flap hinges, but no lag hinges or cyclic pitch) and was fabricated with 
very high torsional rigidity (*>9^^* 22 p). A description of this semi-articu- 
lated model and its performance results are presented in Ref. 10. Selected 
performance results for this semi-articulated rotor are included in these 
figures to show what effect, if any, might result from the elastic character- 
istics peculiar to bearingless rotors. The semi-articTilated rotor resTilts 
presented in these figures are cross-plotted values to obtain results at the 
three advance ratios of .25, .35, and .45, and for tip-path angles of attack. 
This data maniptxLation was required to approximate more closely the bearingless 
rotor model conditions wherein, because of the deliberate ntilling of hub 
moments, the shaft angle of attack coincided with the tip-path plane angle of 
attack. 

Figures 23, 24, and 25 show that the cross plotted resixLts for the semi- 
articulated model rotor generally agree quite well with the bearingless rotor 
results. Some descrepancies in collective pitch could arise due to the 
accuracy with which collective pitch could be set and with which the tip-path 
plane can be measured in the case of the flapping rotor results. The pinned- 
pinned torque tube configuration performance results of these three figures 
can be used, in a like manner as those for hovering performance, as a basis 
for interpreting the cantilevered torque tube configuration resTilts. 

Figures 26 throu^ 29 show the forward flight performance results for 
the cantilevered torque t\ibe configured model at an advance ratio of 0.25 
for each of the four selected push-rod attachment locations. Comparison of 
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these figures with the corresponding results for the pinned-pinned torque 
tube configuration (Fig. 23) shows a reinforcement and extension of the 
hovering performance results: that a significant amount of control angle 

washout occurs with the cantilevered torque tube configurations. Furthermore, 
also in repetition of the hovering performance results, the two trailing 
edge attachment configurations have more washout at the higher collective 
singles than do those configurations with leading edge attachment. Figures 
30 and 31 present similar cantilevered torque tube forward flight performance 
results for an advance ratio of .35> for two of the selected push-rod 
attachment points. Figures 32 and 33 then present corresponding results for 
an advance ratio of .47. In all these figures the same general findings as 
discussed above are reinforced. 

Vibratory Stresses - Figures 34, 35 > and 36 present l/2 FTP values of in- 
board flatwise bending stresses, inboard edgewise bending stresses and outboard 
flatwise bending moments, respectively, for the model with the pinned-pinned 
torque tiibe. Note that these metric system stresses can be converted to the 
English system using Table I. General observations to be made from these data 
are (l) that the l/2 PTP bending moments and stresses are monotonic increasing 
functions of advance ratio, (2) that the flatwise bending stresses and moments 
are monotonic increasing factions of blade loading, C^/ct^ (3) that edgewise 
bending stresses not only Increase faster with blade loading than do flatwise 
stresses, but appear to achieve minimimi values at finite values of blade 
loading, and (4) that the stresses are relatively insensitive to the means by 
which blade loading is achieved (collective angle vs. shaft angle of attack). 
Figures 37> 38 and 39 present corresponding 1/2 PTP bending stress and moment 
results for the model rotor configured with the cantilevered torque tube and 
with one of the trailing edge push-rod attachment locations. Similarly, 

Figs. 40, 4l, and 42 present l/2 PTP bending stress and moment results for 
one of the leading edge push-rod attachment locations. Observations which can 
be made from these figures, when compared with corresponding pinned-pinned 
torque tube results and with each other are: (l) that the inboard vibratory 

flatwise stresses for the cantilevered torque tube configurations showed 
stronger dependency upon blade loading than did those for the pinned-pinned 
torque tube, (2) that the cantilevered torque tube cases showed considerable 
0.5 P oscillatory behavior at the higher advance ratio (p. = 0.47), (3) that, 
for the cantilevered torque tube configuration with trailing edge push-rod 
attachment, the vibratory stresses were sensitive not only to C^/a, but in- 
dependently to collective angle as well, (4) that, at the low blade loadings, 
inboard stresses for the cantilevered torque tube configurations were less 
than those for the pinned-pinned torque tube configuration, but that the re- 
verse was true for the outboard flatwise bending moments, and (5) that the 
edgewise stresses for both torque tube configurations showed close qualitative 
agreement with variations in blade loading. An important characteristic of 


22 



the cantilevered torque tube configurations is that the (vibratory) push-rod 
load directly contributes to the total vibratory bending stresses in the 
flexbeam. Thus, it is not unreasonable to find that the vibratory bending 
stresses should be independently affected by blade control angle. The source 
of the 0.5 P oscillations is unknown, but the results of the correlation study 
suggest that the rotor may have been in a once every other revolution re- 
treating blade stall condition, perhaps precipitated by the inherently soft 
torsion system for cantilevered torque tube conf igiu'ations . 

Control Angle Derivatives - For selected trimmed test conditions all rotor 
control angles (A^g, were independently perturbed to obtain 

the variations in performance and stresses due to variations in these control 
angles. Table VII summarizes the performance perturbation results in the form 
of partial derivatives of the measvtred load coefficients (divided by solidity) 
with respect to each of the control angles. These results again demonstrate 
the loss of control power experienced with the cantilevered torque tube con- 
figurations relative to the pinned-pinned torque tube configuration. Also 
shown in this table is that in this advance ratio range the lift coefficient/ 
solidity derivatives are monotonic increasing functions of advance ratio. 

This table further corroborates Fig. 22 in showing that the cantilevered 
torque tube configviration with leading edge push-rod attachment has signifi- 
cantly more control power than does that with trailing edge attachment. This 
is consistent with the fact that the geometric pitch-flap coupling for the 
leading edge attachment configiiration is positive (destabilizing) whereas that 
for the trailing edge attachment configuration is negative (stabilizing). 
Geometric pitch-flap coupling is defined herein as that calculated using pitch 
horn and mode shape geometry and discounting the washout due to flexbeam 
flexibility effects. The difference in pitch-flap coupling values for the two 
cantilevered torque t\ibe configurations is also seen in the difference in 
sign for the lift coefficient derivative with respect to longitudinal cyclic , 

S (Cj/a)/dAig 

The stress pertubation res\ilts are given in Figs. 43 through 47 for 
the same torque t^ube configurations and conditions given in Table VTI. These 
figures each present the time-histories over one azimuthal revolution of 
inboard flatwise and edgewise bending stresses for the trimmed and plus and 
minus perttirbations of control angles. Most notably these figures show that 
the variations of edgewise bending stresses with perturbations in control 
angles are generally small relative to those for flatwise bending stresses. 

The only case showing significant variations in this stress is the .35 
advance ratio case for the cantilevered torque tube with leading-edge push- 
rod attachment location (Figs. 47a through 47d). Less obvious from these 
figures is the result that the variations in these stresses appear to be pre- 
dominantly first harmonic. Furthermore, some cases (Figs. 44a, 44d, 46d, 
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^7b» 37d) shov significant differences in an^ilitude of variation for the 
plus and minus control angle pertuhation , this vould indicate nonlinear 
functional relationships of these stresses to the control angles. 
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CORRELATION RESULTS AND ANALYSIS VALIDATION 


Time-histories of stresses and accompanying performance data for select- 
ed runs from the experimental model tests were used as correlation data for 
validation of the rotor aeroelastic analysis (computer program GUOO, described 
in Ref. l). Discussion of the details of the test model and of the experi- 
mental wind tunnel test procedures are contained in a previous section. This 
section discusses the correlation and analysis validation study which con- 
sisted of (l) selection of appropriate experimentail test cases for correlation, 
( 2 ) formulation of valid ground rules for meaningful correlation and modifica- 
tion of the analysis to comply with them, and (3) actual validation of the 
analysis by applying it to the experimental test data. The following sub- 
sections discuss these component tasks. 


Selection of Correlation Cases 


A total of twelve (l2) test cases were selected for the purpose of 
program validation. The major criterion for selection was to give a range 
of blade loadings (Cj^/o), advance ratios, and shaft angles representative 
of realistic fiill-scale flight conditions. In addition, cases were specifi- 
cally chosen to substantiate the operation of the program within limited 
paraiaeter variation. Two cases were chosen at substantially the same 
blade loading, advance ratio, and shaft angle but with two and four blades, 
respectively, to establish, if possible, the independence of program correla- 
tion to number of blades. Where possible, both pinned-pinned and cantilevered 
torque tube configured cases were chosen at comparable load conditions, to 
validate program operation for each of these optional dynamic configurations. 
Within the cantilevered torque tube configured cases, both leading- and 
trailing-edge push-rod attachment location cases were selected. Finally, 
since most experimental cases were run with the hub moments (i.e., IP flatwise 
stresses) very close to zero, two cases were selected from a run with 
perturbations in control angles to show the isolated effect of substantial 
hub moment. All of the twelve correlation cases used for the program 
validation are sinmnarized in Table VIII. 

It is to be noted in this table that the pinned-pinned torque tube had 
essentially zero pitch- flat /edge coupling for all modes, whereas the 
cantilevered torque tube configuration demonstrated significant (geometric) 
pitch-flat coupling. Because of the existence of a "wobble mode" (rigid body 
torsion coupled with flatwise bending through the flexbeam bending flexibility) 
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for this configuration, actual pitch- flat coupling is dynamically phase lagged 
and modulated from the geometric value (see Ref. l). The geometric shaft 
angles listed must be interpreted and used with reference to the zero lift 
hub angles of attack; for advance ratios of . 25 , *35j and .UT the zero lift 
hub angles-of-attack were, repsectively : -0.3, -0.7 and -1.0 degree. Finally, 

the control angles indicated in this table are those as best estimated 
from the considerations discussed in the following subsection. 


Groiond Rules for Establishing Correlation 


In establishing the proper basis for correlation, it became apparent 
that the following factors had to be considered: 


1. Using the data of Table VI, it can be shown that, for the outer 
25 percent span of the blades , the various flight conditions tested 
defined a Reynolds number range of 80,i+00 to 422,000' and a correspond- 
ing Mach number range of .08l to .428. For this Reynolds number range 
and considering the clean airfoil construction used on the model 
blades, boundary layer transition can be expected to strongly affect 
the airfoil characteristics. Indeed, Ref. 10 shows, for a rotor 
operating in this test condition - Reynolds number range, that: 

a) the advancing blade sections are exclusively laminar, b) the low 
Reynolds number retreating blade sections are predominantly turbulent, 
and c) the extent of transition is a strong monotonically increasing 
function of blade loading. 

2. Reliable NACA 0012 airfoil data in the operating Reynolds number 
range for these model tests are almost nonexistent. Extrapolation of 
existing (turbulent boundary layer) MCA 0012 airfoil data into the 
low Reynolds number range thus becomes more of an expeditious procedure 
than an accurate estimate. 

3. Because the bearingless rotor is a hingeless rotor, substantial IP 
stresses can easily be generated from IP blade angles-of-attack. 
Furthermore, these IP stresses are generally capable of masking and/or 
diminishing the higher harmonic content of the stresses. 

4. Several sources of error in the measurement of the IP blade angles- 
of-attack exist: 
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a. Blade control angles were measured primarily at the nonrotating 
portions of the swash-plate by means of potentiometers. These 
measurements Inherently lack the ability to account for torsional 
flexibilities in the control angle loeid path above the swash plate. 
Potential sources of such flexibility include the swash plate itself 
and, for the pinned-pinned torque tube configuration the bushing for 
the inboard retention bolt and the Thomas coupling at the outboard 
end. While these latter flexibilities were lacking for the canti- 
levered torque tube, the bending flexibility of the spar provided 

a mechanism for considerable "washout" of the Inputted control angle 
(see Ref. 1 , for a discussion of this flexibility). This represent- 
ed a significant problem in determining the actual blade angle under 
rotating conditions for the cantilevered torque tube configuration. 

b. Reference signals from the potentiometers for the zero values of 

control smgles could at best, be obtained only indirectly: at zero 

rotor speed using an inclinometer or, at nominal rotor speed, using 
thrust measurements. 

c. The effective rotor (hub) angle-of -attack could only be estimated 
from the measured geometric angle-of -attack and the "zero thrust" 
angle-of -attack, i.e., the geometric angle needed to null the thrust. 
This latter quantity essentially represents the vertical deflection 
of the trimmed flow caused by model blockage. At nonzero thrusts, 
however, both blockage effects and timnel-floor reflection effects 
would be expected to increase. These two effects are, furthermore, 
extremely difficult to measure or estimate and no attempt was made 

to do ‘so. 

d. It is likely that significant IP components of inflow exist 

within the test section. Beyond any possible txmnel swirl (existing 
in the absence of the test model) a mechanism exists for asymmetric 
flow conditions: the motor torque which is imparted to the air by 

the rotor creates, in effect, a differential lateral blockage. That 
is, the equivalent drag center for the rotor model is laterally dis- 
placed from the tunnel plane of symmetry. 

The above four items are all probably present to some extent, some greater than 
others. Their uncertainties are attributable to the e:qpense and/or a lack of 
practical method of reliably and accurately measuring these effects. However, 
these effects are not considered InfluentisOL in achieving the primary objective 
of this study which is to evaluate the overall feasibility of the Composite 
Bearingless Rotor. 
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5. Measurement of rotor thrust is relatively reliable and accurate; 
repetitive reference signals for zero thrust were easily obtained re- 
sulting in good calibrations. 

6. The effects of blade pitch angle and structiural twist of the flex- 
beam (resulting from pitch angle) on edgewise loading and edgewise 
coupling with flatwise should be significant based upon the order of 
magnitude of terms in the coupled edgewise bending response equation. 
Hence, the actual blade pitch angles should be as closely duplicated 
in the analysis as possible. 

A consideration of these factors led to the following selected procedures 
for correlation: 

1. The airfoil characteristics to be used as input to the aeroelastic 
analysis were estimated from the results of numerous previous aero- 
dsmamic studies made of rotors tested in this wind tunnel. Zero varia- 
tion in Mach number (and, equivalently, Reynolds number) was assumed 
and a typical angle of attack variation in c^g, c, , c^ /, was assumed. 
Based upon known stall boundaries for similar rotors tested in this 
facility (see Ref. 10) a lift stall airfoil section angle of attack of 
11 degrees was assumed, with a gradual stall behavior above this angle. 
Drag coefficient variation with angle-of -attack was selected to best 
simulate various measured rotor drag and torque measiirements. Moment 
coefficient was selected from low Reynolds number two-dimensional data 
available at UTRC. These assumed airfoil characteristics are plotted 
in Fig. 48 . At best, they are a coarse approximation and represent the 
most probable impediment to accurate correlation; this is especially 
true for the drag coefficient upon which the edgewise stress is so 
dependent . 

2. Meaningful stress correlations could be best made if the analysis 

were first run in a trimming mode wherein the calculated thrust and 
first harmonics of flatwise stresses were appropriately forced to 
agree with the measured values (actual measured thrust and null values 
for first harmonics of inboard flatwise stress). The aeroelastic analysis 
was thereby modified: (1) to trim on hub pitching and rolling moments 

as determined by first harmonics of blade root flatwise stresses, and 
(2) to effect the trim calculation by varying the zeroth and the first 
cosine and sine harmonic components of a generalized Glauert (momentum) 
inflow model (see Ref. 1) keeping control angles fixed. 

3. For the pinned-pinned torque tube configuration, flexbeam torsion 
strain gage data were used to determine control angles. For the 
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cantilevered torque tube configviratlon, however, wherein substantial 
torsion flexibility and pitch-flat bending coupling effects obscure 
the torsion strain gage signals, the swash-plate potentiometer data 
were used, and a rigid-body torsion (wobble) mode was included in 
the analysis. 

4* It is assumed that the primary purpose of the correlation stiidy is 
to validate the aeroelastic analysis in so far as it adequately simu- 
lates the primary characteristics peculiar to the bearingless rotor. 
Therefore, although the analysis has provision for using variable 
(vorticity induced) inflow and unsteady (dynamic stall) aerodynamics, 
it was decided to use a more conventional inflow description (general- 
ized Glauert using momentum considerations) and the usual quasi-static 
airloads formulation. This decision was further Justified in that 
first, the variable inflow would be expected to affect only the hl^er 
harmonics of stress at the cost of significantly increased computation. 
Secondly, since the static airfoil characteristics were known only 
approximately, the inclusion of \msteady data was unwarranted. These 
two simplifying approximations should be kept in mind in the discussion 
of the results of applying the analysis, which follows. 


Application of Test Results to Analysis Validation 

Using the parameters defining the twelve correlation cases given in 
Table VIII together with the considerations discussed in the above subsection, 
the rotor aeroelastic analysis was run to simulate these test conditions. The 
results of these correlation runs are given in Table IX , and Figs. 49 through 
86. For each of the twelve correlation cases, the analysis was first run with- 
out a trim calculation using the control angles listed in Table VIII, and then 
with the modified trim calculation wherein rotor lift and first harmonics of 
flatwise stress at the root were matched to the experimental values. The trim 
was obtained by maintaining the control angles at the experimental values and 
varying the generalized Glauert momentum coefficients, \q ^Is* 

each figure the experimental values, and analytic values with and without trim 
are presented. 

The perfonnance correlation results are given in Table rx. Specifically, 
the rotor lift, drag, and torque coefficients (per solidity) for each of the 
twelve correlation cases are listed for the experimental, analytical (non- 
trimmed), and analytic (trimmed) results, respectively. The conventional 
(uniform) and generalized Glauert momentiim inflow ratios, being aerodynamic 
idealizations, pertain only to the analytic result cases, as appropriate. 
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Comparison of the experimental results with the nontrimmed analytic 
results shows that, indeed, the nontrimmed analytic performance results are 
considerably in disagreement with the experimental performance results. 
However, comparison of the experimental results with the trimmed analytic 
results shows substantial improvement in the drag correlation and moderate 
improvement in torque. The corresponding lift results are, of course, 
essentially identical by virtue of the trim calculation. As was discussed 
in the above experimental model test results section the experimental accuracy 
limits on the lift, drag, and torque coefficient (per solidity) were, res- 
pectively, ±.002, ±.004 and ±.003. Hence, it is reasonable that correlation 
is seen to be best for the higher rotor loaded test cases; with the exception 
of the drag results for case 72.06 the trimmed analytical results correlate to 
within the experimental accuracies given above. 

The stress time-history correlation results are shown in Figs. 49 through 
86. These resxalts consist of time histories of blade flatwise and edgewise 
stresses at the spar root ( r = ,06 R) and blade flatwise moment at a span 
of .42 E, which is somewhat outboard of the torque tube to flexbeam jmcture 
(r = .31 R). For the four cantilevered torque tube cases, results are 
additionally presented for the flexbeam torsional stresses equivalently 
expressed in terms of flexbeam twist deflection. The results for the twelve 
cases defined in Table VIII are presented sequentially in the fig\ares, and 
each of the four Roman numerically designated groups in Table VIII are 
discussed separately below. 

As with the performance correlation results, stress time-histories are 
presented each for the experimental, analytic (nontrimmed) and analytic 
(trimmed) results. 

Group I - Effect of Number of Blades - Since the aeroelastic analysis 
used in this study assmes hub rigidity and momentum inflow, it is inherently 
independent of numbers of blades. Hence the time-histories presented in Figs. 
49 through 54 mainly show the effect of number of blades upon the experimental 
resiilts for the pinned-pinned torque tube configuration. The effect of number 
of blades is essentially aerodynamic and, hence, it is reasonable to expect 
that the results found for one torque tube configuration would apply to both. 
Case 25.03 (Figs. 49 through 51) and case 87. 08 (Figs. 52 through 54) are both 
at approximately the same flight conditions as defined by advance ratio, blade 
thrust loading, and control and hub angles. The lowest advance ratio (n = .25) 
was chosen to accentuate the aerodynamic differences. From comparisons of 
Figs. 49, 50 and 51 respectively, with Figs. 52, 53 and 54 the following obser- 
vations can be made. 

1. The qualitative features of the experimental stresses and moments 


30 


(j)eak-to-peak levels and phase angles) are relatively insensitive to manber of 
blades. 

2. The four-bladed (case 87 . 08) stress time histories show greater har- 
monic content than the two-bladed case time histories. This should reasonably 
follow from the increased variable inflow due to the wake vorticity of the two 
additional blades. 

3 . The nontrimmed analytic root flatwise bending stresses show substan- 
tial lack of correlation with the experimental values. 

4. The effect of trimming is to greatly improve the root flatwise bending 
stress correlation. In general, the phasing and harmonic content (mostly 
second) are well simulated. 

5 . The trimmed analytic root edgewise bending stress restilts show moder- 
ate improvement over those for the untrimmed analytic results. The main im- 
provement appears to be in the phasings of the first harmonic contents; the 
first harmonic amplitudes do not agree well for these cases. 

Group II- Pinned-Pinned Torque Tube Configuration, Flight Condition Variation - 
Time history results for cases 22.03, 22. 04, 21.05 ^nd 22.05 sre shown in 
Figs. 55 through 66 , and show the effects of limited fli^t condition variation 
upon correlation for the pinned-pinned torque tube configuration. The impor- 
tance of this torque tube configuration to the correlation study is that it is 
virtually free of any complicating pitch-flap coupling characteristics in con- 
trast to the albeit mechanically simpler cantilevered torque t\jbe configuration. 

In general, the correlation results for these four pinned-pinned torque 
tube cases show characteristics similar to those observed and discussed above: 
While the untrimmed analytic results correlate very poorly, the effect of 
trimming is to give relatively good correlation. Specific observations which 
can be drawn from the trimmed analytic results of Groups I and II combined 
are as follows: 

1. Root flatwise stresses (Figs. 49, 52, 55, 58, 6 l and 64): The pre- 

dicted second and higher harmonic contents and phasings are generally 
in good agreement with experiment. Correlation is poorest at the low 
advance ratio cases where the observed increased higher harmonic con- 
tent is not simulated; correlation of harmonic content is consistently 
excellent at the two higher advance ratios. The mean or zeroth harmonic 
stresses do not consistently correlate well, however. 
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2. Root edgewise stresses (Figs. 50, 53> 56, 59, 62, and 65): These 

stresses are both observed and calculated to be predominantly IP and the 
predicted values show moderate agreement in amplitude and phase with 
experiment. The predicted values are generally smaller in PTP amplitude 
than, and are phased more in advance of, experimental. These cases re- 
present concurrent variations in both blade loading and advance ratio. 
Coinparison of the results shows that the azimuthal maximum and minimum 
points of the IP component appear to be sensitive to these two parameters. 
While not shown in these correlative experimental resiilts, the effect of 
blade loading change at constant advance ratio was examined using harmonic 
analyses of stresses at other flight conditions. There, it was consis- 
tently seen that at a given advance ratio the phasing of the IP content 

of the edgewise stress was monotonically advanced (i.e., maximum stresses 
occur at an earlier azimuth angle) with increasing blade . loading or 
equivalently, with mean airfoil angle-of-attack. Thus, since the cases 
of Group II are also at monotonically decreasing blade loadings, the trend 
of IP edgewise stress phasing seen in these correlation cases is consis- 
tent. Although the analytic edgewise stresses are, as noted, generally 
phased in advance of the experimental results, the combined blade loading- 
advance ratio trend appears to be reasonably well correlated. More will 
be said on this topic of IP edgewise stress phase angle below in the 
discussion of the next group of correlation cases. 

3. Blade outboard bending moment (Figs. 51, 5^, 57, 60, 63, and 66); 
Despite a substantial variation in advance ratio and blade loading the 
experimental results show the remarkable similarity that the responses 
are predominantly IP, 3P, and 4p, and the peak-to-peak values and general 
azimuthal signature are all invariant. The theoretical results show a 
similar qualitative behavior in that there is no significant PTP variation 
with these parameters either. The analytic PTP values however are con- 
sistently less than the experimental values as a result of the analytic 
results being predominantly only IP, 2P, and 3P* 

Group III - Effect of Hiib Moments - Because of the proximity of the first edge- 
wise mode frequency to IP (toy^ = 1*35), the low inherent damping of this mode 
and the potential for significant elastic coupling of edgewise bending with 
flatwise bending due to twist, it is reasonable to allow that blade edgewise 
IP stresses mi^t be a function of IP flatwise bending, or equivalently rotor 
hub moment. Since the aeroelastic analysis predicts zero elastic coupling of 
flatwise and edgewise bending due to twist a further test of the analysis 
is to compare analytically predicted edgewise stress responses with experimen- 
tal. ones to detect this coupling, if present. In order to isolate this effect 
correlation was attempted with cases 11.05 and 11. 06 which differ by 1.3 
degrees of B, as shown in Table VIII. Case 11.05 is essentially a zero hub 
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moment (zero IP flatwise *bending) case whereas significant IP flatwise bending 
is present in case 11. 06 as is shown by contrasting Figs. 67 and 69. Figures 
68 and 70' show the corresponding edgewise stresses for these two conditions. 

It can readily be seen that, although the analytic results show a substantially 
advanced phase angle from the experimental results, the FTP edgewise stresses 
and the neglibility of IP flatwise bending upon edgewise bending are well 
correlated. As was discusssed in the Group II results the phasing of the IP 
edgewise stresses appears to be a strong function of some combination of 
blade loading, or equivalently, mean blade section angle of attack and of 
advance ratio. Indeed, cases 11.05 and 11. 06 being low disk loading cases, 
show a relatively retarded IP phase angle and are thus consistent with this 
general trend discussed earlier. 

Group IV - Cantilevered Torque Tube Configuration, Flight Condition Variation^ 
Unlike the pinned-pinned torque tube, the cantilevered torque tube intro- 
duces a significant degree of elastic coupling between flatwise bending and 
rigid body feathering. The primary resulting coupled mode is referred herein 
as the "wobble" mode. Direct consequences of this coupling are: a) control 

angle "washout" wherein swash-plate input angles, considerably larger than those 
for a pinned-pinned configuration, are required to achieve comparable perfor- 
mance and/or trim, b) general inability to eliminate pitch-flat coupling for 
more than one bending mode at a time, c) resulting pitch-flat couplings mani- 
festing dynamic amplifications and phase lags due to the impedance character- 
istics of the wobble mode, and d) a relatively "soft" torsion system with the 
potential for coupling with all the major lower frequency modes. 

A mathematical requirement for the successful modeling of this elastic 
coupling is the ability €0 calculate the force-deflection characteristics of 
the flexbeam- torque tube assembly at the push-rod attachment point. That is, 
the out-of-plane deflection of this attachment point accruing from blade flex- 
beam bending flexibility and due to a parallel directed unit load at that point 
must be accurately calculated. Reference 1 describes some of the details in- 
volved in this calculation, which is an internal portion of the aeroelastic 
analysis and uses transfer matrix techniques. This calculation was satisfac- 
torily validated experimentally by measuring the flexibility of a statically 
mounted blade spar (flexbeam) loaded in tension to simulate the centrifugal 
load. For the nominal condition, a flexibility of mm/N was calculated 

versus a measured value of *hlk mm/N for the blade spar with an equivalent 
static torsion load. The use of this flexibility for modeling the elastic 
restraint and "dynamic" pitch-flat coupling of the wobble-mode is also des- 
cribed in Ref. 1. 

The time-histories, shown in Figs. 71 through 86, present the correlative 
results for the cases of Group IV. Specifically, they show comparisions between 
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experiment and theory for limited variations in flight ^conditions and 
pitch-flat coupling, as determined by push-rod/torque tube attachment point 
location. In addition to those stress and moment quantities similarly shown 
for the pinned-pinned torque tube configuration, flexbeam torsion deflec- 
tion results are also presented. The torsional rigidity of the pinned-pinned 
torque tube configuration discussed above in Groups I, II, and III was rela- 
tively high (i «0 = 6 . op). Hence, the flexbeam was effectively vincoupled 

in torsion so that its torsion stress, or equivalently, elastic twist was 
substantially determined by the swash-plate input angles. Indeed, the time- 
histories of the flexbeam torsion deflections for the pinned-pinned torque 
tube cases were regularly sinusoidal and consequently used to establish the 
effective control input angles. Because of their regularity, formal presen- 
tation of these time-histories would have been superfluous and was consequently 
omitted. In the case of the cantilevered torque tube, however, the torsion 
system was relatively soft, owing to the inclusion of flatwise bending flex- 
ibility in the effective torsion restraint. The time-histories of flexbeam 
elastic twist, in this case, became multi-harmonic; superimposed on the l/rev 
character due to swash-plate inputs were the multi -harmonic responses of the 
wobble-mode. Thus, the flexbeam torsion traces ceased to be a source of input 
calibration and, moreover, became an additional correlative quantity. Specific 
observations to be drawn from the trimmed analytic results of Group IV are as 
follows : 


1. Blade root flatwise bending stress (Figs. 7 I 5 75? 79> and 83 ): The 

predicted results follow the same trends observed for the pinned-pinned torque 
tube. In general, the low advance ratio (p = .25) resiilts underpredict the 3P 
(and higher) resiilts, but the higher advance ratio (p = 0.47) results are in 
excellent agreement. The mean value is also well correlated but not consistent. 

2. Blade root edgewise bending stress (Figs. 72, 76, 80 and 84): Again 

the predicted results follow the same trends observed for the pinned-pinned 
torque tube. The predominantly IP FTP amplitude is underpredicted and phased 
in advance of the experiment. 

3 . Blade outboard bending moment (Figs. 73? 77, 8l and 85 ): Whereas the 

lower advance ratio results are seen to parallel the results for the pinned- 
pinned torque tube, the predicted harmonic results for the higher advance 
ratio case are in very good agreement with the experiment. 

4. Flexbeam torsion deflection (Figs. 74, 78, 82 and 86): The predicted 

results are predominantly IP and correlate the IP in amplitude and phase 
excellently. It should be noted that these time-histories each represent the 
sum of collective and cyclic swash-plate inputs together with the substantial 
wash-out due to mean and IP wobble mode responses, as reference to Table VIII 
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shows. The higher harmonics (mainly 2P, 3P, and 4p) seen in the experimental 
results are not well simulated; this is consistent with the underprediction 
of both 3P root flatwise bending and IP root edgewise bending whose product 
would thus constitute an undeipredicted 4p torsion response. 

5. Subharmonics of case 73.08 (Figs. 83 thru 86 ): These figures each 
show two experimental time histories for the plotted azimuthal period of 360 
deg. : these multiple traces each represent the two halves of time histories 
which repeat in a regular fashion with an azimuthal period of 720 deg. Thus, 
the responses are manifesting (stable) O. 5 P subharmonic characteristics. 
Examination of the experimental portion of these figures shows that: 

a) The experimental 0.5P subharmonic content is most dominant in the 
blade root flatwise and edgewise bending responses, 

b) The experimental flatwise bending stresses (Fig. 83 ) show the 
greatest degree of dissimilarity between revs on the retreating to 
aft portion of the disk. The dissimilarity is one primarily in level 
and secondarily of phase, since the maximum and minimum points occur 
at approximately the same azimuth angles. 

c) The dissimilarity of edgewise bending stresses (Fig. 84) between 
revs is mainly one of phasing; the transition from maximum stress to 
minimum is either lagged or advanced between revolutions, again on the 
retreating to aft portions of the disk. 

d) On the aft, advancing to leading portions of the disk the experimental 
flatwise bending and torsion responses both show a dissimilarity between 
revs in 2P amplitude with the phasing virtually unchanged. 

The above results can be interpreted to a certain extent by noting that 
in Fig. 83 , on the retreating side for revolution no. 1, the flatwise bending 
response is unable to sustain the increased 2P response amplitude shown on 
the advancing side. Moreover, at an azimuth of approximately 225 degrees, the 
rev. no. 1 appears to demonstrate a stalled characteristic. Furthermore, 
examination of Figs. 84 and 85j respectively, show for the rev. no. 1 portions 
of the responses, that at the forward portion of the disk the blade concur- 
rently has an increased aft lag rate and an increased pitch angle, both con- 
ductive to stall. The analytic predictions on the other hand show a similar 
stalled characteristic. Indeed, although not shown, examination of the details 
of blade angle of attack distribution for this analytic case show the blade to 
be mildly (quasi-statically) stalled over the azimuthal sector defined from 
230 degrees to 330 degrees. Thus, it would appear probable that a principal 
source of the experimental 0.5F subharmonic response is a subharmonic stalling 
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of the hlade on the retreating side. It is likely that this stalling is 
aggravated by the subharmonic increase in 2P torsional responses which are 
a probable result first of the higher coupling of torsion with the flatwise 
bending responses, which themselves are mainly 2P, and second, the greater 
susceptibility of blade torsion to (relatively) low frequency excitation 
due to the increased torsional flexibility of this torque tube configuration. 

Correlation Study Evaluation & Comments 

Based upon the details of the correlation cases discussed above and 
recalling both the experimental limitations and the assumed aerodynamic 
simplifications, it is generally concluded that the aeroelastic analysis 
predicts blade bending and torsion responses consistent with rotor 
performance with reasonable accuracy. It is fiorthermore concluded that 
the analysis is suitable for the aeroelastic analysis of the full scale 
design. Other more specific comments to be made concerning analysis 
correlation are as follows: 

1. It is likely that the analysis was run with inaccurately defined 
flight -condition parameters for several of the correlation cases in that the 
analysis of some of the correlation cases ( 87 . 08, 22. OU, 21.05, and 73.08) 
achieved exceptionally good correlation for all stress and moment variables. 

2. Correlation tended to improve with advance ratio. This trend is 
most likely due to the known reduction of variable inflow effects with speed. 

3. Best correlation was achieved with Inboard flatwise stress. The most 
probable reason for this is that this variable is so strongly a function of 
lifting airloads which are the most linear and can thus most accurately be 
approximated. 

Actual one-to-one correlation of edgewise stress both in magnitude 
and phase was only moderate. However, the trends of edgewise stress with 
parametric variation were duplicated with reasonable accuracy. Since the 
edgewise stresses were predominantly one per rev it is most likely that the 
edgewise stresses were governed mainly by aerodynamic drag forces. It is 
therefore likely that the only moderate correlation of edgewise stresses is 
due to the coarseness of the assumed model airfoil drag characteristics. 

5 . Correlation of the predominantly one/rev flexbeam torsion stresses 
for the cantilevered torque tube was excellent. Higher harmonics of this 
torsion stress were consistently underpredicted, however. 
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6. The 0.5P subharmonic responses of case 73 •OS were not successfully 
simulated* One probable reason for this discrepancy is that the model rotor 
appeared to be stalling every other revolution and was therefore a strongly 
aerodynamic phenomenon requiring accTorate airfoil drag data which wasn’t 
available . 

7» The analysis accurately predicted a minimal influence of IP flatwise 
bending on edgewise stresses. 

8. The analysis correlated well in a passive sense in that it did not 
predict \mstable motion for any case which was observed to be stable. 



FULL SCALE DESIGN 


The primary objective of conducting a preliminary full-scale design as 
part of this feasibility study were 1) to enable an experienced rotor designer 
to examine the practical aspects of the CBR from the viewpoint of fabrication 
and total system requirements, and 2) to provide realistic blade mass and stiff- 
ness characteristics for use in the aeroelastic analysis discussed in a later 
section of this report. The specific rotor size was not of particular impor- 
tance. It was believed, however, that a relatively large rotor would be of 
greater significance since there has been a general concern over the ability 
to scale up relatively small hingeless rotors (diameters less than 12m (40 ft)) 
for large helicopter applications. An S-6l size aircraft was selected having 
a 4-bladed rotor with a diameter of l8.9 m (62 ft), chord of .58 m (1.9 ft), 
and a design tip speed of I98 m/sec (65O ft/sec). Initially, blade frequencies 
typical of conventional hingeless rotors were selected, i.e., first mode flat- 
wise frequency slightly greater than one-per-rev and first edgewise frequency 
between one-per-rev and two-per-rev. It was felt that such a design would be 
relatively free of instabilities and the aeroelastic analysis could focus ex- 
cl'uslvely on the dynamic characteristics associated with the CBR and not be 
influenced by unrelated d3mamic phenomena of hingeless rotors. Two designs 
were considered, one having a cantilevered control torque tube, and the other 
having a control torque tube with a "snubber" support at the root end to elim- 
inate pitch coupling associated with flapwlse motions of the torque tube. The 
decision to Include the snubber design was influenced primarily by an early 
wind tunnel investigation in which adverse coupling was experienced under cer- 
tain low tip speed operating conditions. Although the test was exploratory in 
nature and the coupling arising from an unrealistic off-design operating con- 
dition, it demonstrated the potential for undesirable coupling in the canti- 
lever torque tube design which should be thoroughly evaluated before proceeding 
to full-scale development. It was thus decided that both torque tube concepts 
be configured so that preliminary design information would be available for 
each. 


A comprehensive aeroelastic analysis of the full-scale design is presented 
in a later section, and in this analysis only the "snubber" design is studied. 
The extent of this study was such that only one design could be so investigated, 
and it was decided that the first generation CBR to reach flight status would 
probably be of the snubber design — the minimimi risk approach. However, both 
designs were evaluated in the model wind tunnel program. 

Several composite materials were available that could possibly satisfy the 
stiffness requirements of the CBR. Most notable were fiberglass, boron/epoxy, 
graphite/epoxy . From the standpoint of cost, fiberglass was the. obvious 
choice; however, due to its lower bending modulus, ein increased flexbeam 


38 



cross-section would be required to achieve the desired frequencies. This wotild 
also necessitate an increase in flexbeam length since its torsional stiffness 
would go up due to the larger cross-section, fhe shear modulus of fiberglass 
is approximately equivalent to that of graphite /epo^Qr so that no benefit in re- 
duced torsional stiffness can be realized in using fiberglass with low bending 
modulus. Graphite/epoxy was selected over the boron/epoxy primarily due to its 
lower cost and greater machinability. It was also favored because of the ex- 
perience gained in its use during the initial studies and experiments conducted 
in developing the CBR concept. 

The design features of both the cantilever and the snubber torque tube 
configurations are described in the following paragraphs. 


Cantilever Torque Tube Design 

The general layout of the cantilever torque tube design is shown in Pig. 
87 . The initial requirements that set the basic spar (dimensions were a first 
flatwise mode frequency between 1.05 P and 1.10 P to ensure good handling 
qualities, a first mode edgewise frequency greater than 1.25 P to provide ade- 
quate separation from the primary edgewise excitation frequency of 1 P, and a 
flexbeam torsional stiffness sufficiently low so as not to exceed the push-rod 
design loads or the control actuator limits. The flatwise and e(igewise fre- 
quencies are dete3miined primarily by the flexbeam thickness and chord dimen- 
sions. The torsional stiffness is determined primarily by the flexbeam 
thickness and length. The .flexbeam that satisfied these requirements was 2.44 
m (96 in.) in length, tapered in chord from 0.305 m (12 in.) at the hub 
(station 15) to 0.203 m (8 in.) at Station 111, and tapered in thickness from 
3.05 < 2 m (1.2 in.) to 2.54 cm (1.0 in.) over this same span. The flexbeam 
could either be Joined to the blade at Station 111 by means of a conventional 
bolted Joint or carried through to the tip forming the primary structure of the 
entire blade. The decision to carry through to the tip or have a Joint at 
Station 111 depends upon composite materials manufactTuring facilities, trans- 
portability requirements, and blade folding requirements. The design selected 
for this stu(3y was of the continuous spar type where the spar carried fr(xn the 
tip of one blade through to the tip of the opposing blade. The chord (dimension 
of the outboard spar was maintained at the value at Station 111 (0.203 m)j 
however, the reduced flatwise stiffness contribution reqaaired of the outboard 
spar allowed the thickness to be tapered from 2.54 cm (1.0 in.) to 0.64 cm 

( 0.25 in.). 

The blade cover is comprised of graphite/epoxy with + 45 deg fiber align- 
ment to provide maximum torsional stiffhess. The skin thickness was set at 
0.102 cm (0.040 in.) which was sufficient to provide a fixed root torsional 
frequency of greater than 5 P* The lea(3ing edge is protected with a steel/ 
nickel plated abrasion strip and extends from the outboard end of the torque 
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tube to the blade tip cap. A small abrasion strip sQ.so covers the leading edge 
of the tip- cap. Nomex honeycomb is used to fill the cavity between the spar 
and cover and structurally bonded to both to form a blade which is relatively 
rigid in torsion. The leading edge incorporates a molded counterweight con- 
sisting of lead shot and adhesive to mass balance the blade at the quarter 
chord. 


Ihe torque tube is cantilevered from Station m inboard to station 19 . 

The torque requires not only high torsional rigidity to transmit control input 
to the flexbeam but also requires high flatwise bending stiffness to prevent 
large bending deflections caused by the pushrod loads. To accomplish these 
requirements both + 45 deg plies and 0 deg plies are laminated resulting in a 
maximum skin thickness of 1.52 cm (.6 in.) at the Juncture tapering to 0.25 cm at 
Station 19 . The thickness also tapers to 0.25 cm at the leading and trailing 
edges of the torque tube. In order to accommodate a total pitch travel of -10 
deg to +35 deg and to allow for the bending deflection of the torque tube 
under the maximum push rod load a total thickness of 43.2 cm (17 in.) was 
required. 


The flexbeam design was initiated by prescribing that the flexbeam 
cyclic shear stress not exceed the design allowables under maximum cyclic 
pitch condition. The alJLowable stress was set at 24 MW/m^ (3500 psi). This 
level was established from the fatigue tests reported in a separate section of 
this report, wherein an average fatigue limit of 48 MN/m^ (7000 psi) (defined 
as a 10 percent reduction in modulus after 10^ cycles) was measured. Applying 
a safety factor of 2 gives the design limit of 24 Mw/m^. The cyclic pitch 
reqtiired at the trimmed high speed condition of 170 kn was approximately 10 
deg. A conservative value of 15 deg was selected to include maneuvers. Using 
the midflexbeam station as representative, the stress is calculated from the 
following equation 



where t and I are the flexbeam thickness and length, respectively. Solving 
for flexbeam length, 

_ ^ 0.0279 X 620. 7 . X 0.262 = 1,89 m (7I+ in.) 

fs ■ 

Further conservatism was added by extending the length to 2.44 m (96 in.) to 
cover possible stress concentrations and variations in G. For this flexbeam 
length the calculated maximum cyclic stress at the midpoint is I8.6 MN/m^ 
(2700 psi), well below the fatigue limit. 
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To minimize the steady torsional stress in the flexbeam a b\jilt-in twist 
of 13.5 deg was used. This represents sin aversige collective pitch estimated 
for a typical helicopter mission. The built-in twist is incorporated during 
the material layup of the flexbeam where the actual mold in which the beam is 
cured is twisted to the required angle prior to layup. The aerodynamic twist rate 
of the outboard portion of the blade is -4 deg, and the airfoil section WACA 
0012. If desired, the twist and airfoil section could be easily modified in 
the mold designs for the spsir and skin; however, for the purposes of this 
feasibility study conventional airfoil and twist were used. 

A preliminary aeroelastic analysis, performed after the preliminary design 
study, revealed lowly damped coupling between several bending modes and the 
rigid body torsion mode (wobble mode acting about pushrod support point) at a 
high forward speed condition. This condition was improved with the addition 
of an 11.35 kg (25 Il3) tip weight. 

The blade mass and stiffness characteristics are given in Figs. 88 thru 90 
and torque tube and flexbeam information is given in. S’igs. 91 thru 93 • A 
summary of the blade uncoupled natural frequencies is presented in Figs. 9^ 
and 95- It should be noted that these figures contain results for both sub- 
critical and supercritical designs. These represent the two configurations 
studied in the aeroelastic analysis and refer to the frequency placement of 
the first edgewise mode. A subcritical designation represents a first edgewise 
frequency of greater than the rotational frequency and is the design described 
in this section. The supercritical represents an edgewise frequency less them 
rotational frequency. The decision to examine the superciriticsil case was made 
following the conpletion of the full-scale design effort described in this 
section. Further discussion of the supercritical case is presented in the Aero- 
elastic Analysis Section. 


Snubber Torque Tube Design 

The layout for the snubber torque tube design is shown in Fig. 96 . The 
general design feattires of this configuration are identically equivalent to 
those of the cantilever design except for the torque tube characteristics. In 
this design an additional support point is provided for the torque tube. In 
addition to the juncture at the outboard end of the -flexbeam and the attach- 
ment point of the push rod, an elastomeric snubber is incorporated at the 
inboard end of the torque tube and bonded to it and the flexbeam . With such 
a three point support essentially all pitch coupling due to flapping is elim- 
inated. In this design the torque tube is not required to transmit bending 
moments, only torsional moments, and thus only ±45 deg fibers are used. The 
skin thickness is constant at 0.25 cm (0.1 in.). Also, to achieve the condi- 
tion of low bending moments in the torque tube, the total thickness of the 
torque tube at the jimcture is reduced to that of the blade itself and the 


4l 


fiber angles are increased beyond db deg to aiffect a reduction in bending 
modulus at this point. iHie resulting structure iremains relatively high in 
torsional stiffness but reduced significantly in flatwise bending stiffness. 

Ihe total torg,ue tube thickness at Station 19 to accommodate the 49 deg pitch 
requirement is 38.6 cm (15.2 in.). The structural and frequency information 
presented in Figs. 88 to 94 also apply to the snubber torque design. 

me total blade weight amounted to 120 kg (265 Ih) and the hub plates 
145 kg (320 lb), m comparing GBR weight with a conventional articulated rotor 
of the same diameter and solidity, an additional 45.4 kg (100 lb) are added to 
account for a somewhat lau^ger shsift to accommodate the larger hub moment capa- 
bility. mis results in a total rotor weight plus added shaft weight of 67I kg 
(l480 lb) which is approximately 23 percent less than the articulated rotor. 
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AEROELASTIC ANALYSIS OF 
FULL-SCALE DESIGN 


\ Selection of Configuration Details and Variants 

The specific component structijral elements of the full-scale design "were 
sized by the fundamental requirements of strength (torque tube and flexbeam), 
and by satisfying basic vibrational characteristics (i.e., absence of integral 
harmonic resonances). Variants of the basic design, however satisfactory they 
might be on the basis of these fundamental requirements, must show satisfactory 
performance in a variety of more comprehensive aeroelastic categories. In 
addition, even within the confines of the basic design, as defined in the 
previous section, considerable latitude exists in the proper selection of 
various details which could profoundly affect the aeroelastic characteristics. 
Consequently, a proper design analysis should include an investigation of the 
effects of parametric variation of these details. To satisfy these require- 
ments the general aeroelastic analysis described in Reference 1 and validated 
in the correlation section of this report was employed. This section presents 
the results of applying this analysis in its various specific categories to 
variants of the full-scale rotor design. More specifically, two basic nominal 
blade designs were selected: one, a subcritical configuration wherein the 

uncoupled first edgewise modal frequency was set near to 1.3P and the other, a 
supercritical case wherein the edgewise frequency was set near to 0.8P. The 
former value was selected for subcritical operation as it provided adequate 
separation from IP and 2P and generally represented the state-of-the-art. The 
latter value was selected as being a conservative value from an air resonance 
standpoint based upon the results of References 11 and 12. 

As was discussed in the previous section, although two different torque 
tube configurations were designed, only the snubber design was used herein. 

The nominal torque tube properties for the snubber design are shown in Figures 
91 and 92 and compared with those for the cantilevered design in Figure 91* 
Initial aeroelastic calculations showed that, for the supercritical design, 
the nominal torque tube stiffnesses raised the effective (coupled) first edge- 
wise modal frequency to 0.88P. Indeed, calculated time history responses 
indicated "excessive” IP edgewise responses for this configuration. Conse- 
quently, the supercritical design was reconfigured with flatwise and edgewise 
torque tube stiffnesses equal to 25 percent of the nominal values as used in 
the subcritical design. This reduction in stiffness lowered the effective 
first edgewise frequency to 0.81P. In addition, whereas the nominal push-rod 
radial attachment point of 21 in. is used for the subcritical design, a value 
of 19 in. is used for the supercritical design. This value was selected to 
eliminate all significant amounts of geometric pitch-flat and pitch-edge 
coupling from the supercritical design. Conventional NACA 0012 airfoil data 
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were used with an incremental drag coefficient of O.CX32 to account for 
operational surface roughness. The built-in precone and prelead angles were / 
assumed to be 3*5 degrees and zero, respectively. In both rotor configure- / 
tions three (3) flatwise, two (2) edgewise, and one (l) torsion (uncoupled) 
modes were assumed for the blade elastic description. A structural damping / 
equivalent critical damping ratio of 0.01^ was assumed for the edgewise modes. 
Except where noted, the full redundant analysis capability of the aeroelastic 
analysis (as described in Reference 1) was utilized. The following subsections 
will describe, in turn, the results of the eigensolution analysis, the choice 
and calculation of the basic trim cases, aerodynamic stability derivatives, 
transient responses due to partial material failure, the effects of unsteady 
aerodynamics, and the results of blade transient stability analyses of both 
the cantilevered torque tube configured test model rotor and the full-scale 
design. 


Eigensolutions 

The eigensolutions of interest that are available are the "vacuum" and 
"nonvacuum" calculations wherein the perturbational airloads are, respectively, 
neglected and included. The "vacuum" eigensolutions are useful in establishing 
coupled modal frequencies and mode shapes due to such coupling effects as blade 
control pitch angle, twist, precone and prelead angles, torque tube stiffnesses, 
section mass and tension center offsets and structural damping. Such eigen- 
solution information is most useful for insuring against integral harmonic 
resonances . The "nonvacuum" eigensolutions are useful for identifying flutter 
and/or aeromechanical instability characteristics in hover, and for providing 
a crude indicator of stability trends in forward flight. The nonvacuum eigen- 
solution considered herein is of the "frozen azimuth" type wherein the peri- 
odicity of the equation coefficients defining the mass, damping and stiffness 
matrices is neglected. Such an analysis is thus incapable of accurately iden- 
tifying stability levels of rotors in forward flight. While recourse could 
have been made to extending the nonvacuum eigensolution to a Floquet theory 
type analysis, such as are alternatively described in References 13, l4, and 15, 
the time-history solution of the fully nonlinear formulation was deemed intrin- 
sically more accurate and relevant to the stability analysis of a trimmed, 
realistically configured rotor blade design. A more detailed description of 
this stability analysis is contained below in a subsequent subsection. The 
remainder of this subsection describes first, the results of using the vacuum 
eigensolution to obtain coupled frequency and mode shape information with 
limited parameter variation and second, the results of using the nonvaccum 
eigensolution for a series of hovering flight conditions. 

Using the vacuum eigensolution, the effects of tip speed, collective 
angle and torque tube stiffness variations on the coupled frequencies and mode 
shapes were determined for the two basic blade designs (subcritical and super- 
critical). For each of these designs the above three parametric variations 
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were made about the nominal conditions defined by a tip speed of I 98 m/sec, 
a collective angle of 12 deg and the snubber torque tube stiffnesses shown in 
Figure 92. 

Coupled (Vacuum) Modal Frequencies. - The effects of the parametric 
variations on the coupled modal frequencies are shown in Figures 97 through 
99 the subcritical design and Figures 100 through 102 for the supercriti- 
cal design. These results were obtained by disregarding all aerodynamic terms 
in the general eigensolution. Figure 97 shows for the subcritical design the 
effect of tip speed variation on the coupling of the primitive "uncoupled” 
modal frequencies (those shown in Figures 9 ^ and 95 )• Not unexpectedly, the 
intersection point of the second flatwise (2F) mode with the first edgewise 
(1E) mode at a tip speed of approximately 75 m/sec shown in Figure 95 becomes 
two diverging branches at that tip speed as shown in Figure 97 - At such a 
point each branch is characterized by comparable amounts of both (2F) and (lE) 
modal content. At the nominal rotor speed the original (lE) branch has become 
predominantly (2F). Note that in all the coupled frequency figures, the pre- 
dominant modal contributors to the coupled modes at the high end of the par , 
eter scale are shown in parentheses, in order of decreasing percent contrib 
tion. Of significance is the result that the first edgewise mode (1 e) cm- 
tributes to such a large degree to both the second and third coupled fla w 
modes. In addition. Figures 98 and 101 both show that for the collect i^/e 
pitch angle range most typical of helicopter operation all modes except j 
third coupled flatwise mode are reasonably well separated from integral h .r- 
monies. Since flatwise modes are relatively well damped and 5P excitations of 
this mode are low, the close proximity of this modal frequency to 5P is not 
considered critical. 

Coupled (Vacuum) Mode Shapes - Due to the variety of coupling mechanisms 
considered in the analysis (pitch angle, twist, coning, structural damping, 
etc.), even in the absence of aerodynamic loads the coupled mode shape compo- 
nents are generally nonuniformly phased. That is, there is no point in the 
oscillation cycle of any vibrating coupled mode wherein all portions of the 
blade concurrently have zero kinetic (or, alternatively potential) energy. 
Hence, coupled mode shape information must include, in addition to the usual 
component deflection distributions, the concurrent component velocity distri- 
butions as well. 

Presented in each of the Figures 103 through I 06 , and for each of the 
six ( 6 ) coupled frequencies typically shown in the preceding figures (97 
through 102 ), are the spanwise distributions of modal deflection and (nondi- 
mensional) velocity in each of three ( 3 ) components of motion for ^ = 

0 deg and 12 deg. The components of motion so depicted are inplane* 
and out -of -plane translation (y^ and z^> respectively), and torsion ( 0 ). 
Variation of the coupled mode snapes with the other parameters ,(tip speed and 
torque tube bending stiffnesses) was found to be qualitatively quite similar 
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to the variations •with collective angle. Results for the subcritical design 
at 0 and 12 deg collective angle are shown in Figures 103 and 104, respec- 
tively, and the results for the supercritical design are shown in Figures 105 
and 106. 

A comparison of the zero collective cases (Figures 103 and 105) shows 
two main dissimilarities: (l) the y^ and velocities associated with the 

respective coupled (IF) and (1E) modes are oppositely phased for the two rotor 
types, (2) the presence of significant pitch-flat coupling for the subcritical 
rotor design and negligible amounts for the supercritical one is reflected in 
the respective amounts (or absence) of torsion motion in the coupled flatwise 
modes. Note that the torsional components of these coupled modes shown in 
Figure 103 clearly demonstrate the locally exaggerated torsional flexing of 
the flexbeam as is required to accommodate the automatic pitch change due to 
flatwise bending. 

The coupling influence of collective pitch angle is shown in comparisons 
of Figure 103 with 104 and/or Figure 105 with 106. This coupling is seen to 
increase the nonuniform phasing of the resultant modes; in general, the com- 
ponent displacement distributions are increasingly dissimilar from their 
respective component velocity distributions . A comparison of Figures 104 and 

106 shows, for both the subcritical and supercritical designs, a high degree 
of blade motion in all components for the cgupled second edgewise bending 
(2E) mode. The significantly larger 9 and 9 content in the coupled mode for 
the subcritical configuration is due mainly to the closer proximity of the 
torsional frequency (10.42P) to the subcritical uncoupled (2E) modal frequency: 
(7.28 p) vs. (4.14p) for the supercritical configuration. 

(Non'vacuum) Coupled Mode Damping Levels - The results obtained from 
limited usage of the nonvacuum eigensolution calculation are shown in Figures 

107 through 110. Each of these figures shows the variations of the critical 
damping ratios for the six (6) coupled mode eigenvalues, with concurrent 
variations in collective pitch angle and (spanwise) variable inflow. The 
inflow (or, equivalently, induced velocity) assumed is conventional (see 
Reference l6) and, at each radial station depends upon the local total pitch 
angle. Figures 107 and 108 present coupled modal damping ratios for the 
subcritical and supercritical configurations, respectively, wherein linear, 
unstalled NACA 0012 section aerodynamic coefficients are assumed. That is, 

C£ is assumed proportional to angle -of -attack using a Mach number dependent 
lift curve slope, c^ is assumed equal to c^ (again with Mach number depen- 
dence) and Cjj. /. is assumed zero. With these simplifications the analysis is 
aerodynamically compatible with other similar existing analyses (References 
17 through 19> for example) and provides a basis for evaluating the effects of 
stall. It is thus not unexpected that the results of the present analysis 
should predict, as do References 17 and l8, flap-lag instabilities and/or 
general reductions in damping at large values of collective pitch angle for 
both subcritical and supercritical rotor configurations. The major 
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instabilities are seen, to involve motions predominately comprised of first 
^flatwise bending and first edgewise bending. Comparison of these two figures 
shows that, for the subcritical configuration, it is the first edgewise (lE) 
eigenvalue 'branch" which goes unstable whereas for the supercritical configu- 
ration, it is the first flatwise (IF) branch. A further comparison shows 
that, for the subcritical configuration, as the instability is approached all 
flat and edge modes lose damping and the branch which goes unstable is already 
lightly damped. In contrast, for the supercritical configuration, the 
- (stable) first edgewise branch increases in damping, but the initially heavily 
damped first flatwise branch precipitously loses substantial damping to go 
unstable. The present findings are qualitatively consistent with the findings 
of Reference 17 which, for subcritical config\irations , predict the lead-lag 
(i.e., first edgewise) eigenvalue branch to be the one to go unstable. Direct 
quantitative comparison of the present results with those of the above cited 
works is difficult at best. Significant structural details (highly nonlinear 
twist, several intercoupling modes, generally , nonuniform section properties, 
etc.) of the realistic blade configuration considered herein are either not 
treated with sufficient rigor or are neglected entirely in these references . 

A more realistic set of nonvacuum eigensolution results is presented in 
Figures 109 and 110 again for the subcritical and supercritical configurations, 
respectively. For these results the complete (nonlinear) conventional MCA 
0012 airfoil data, including both stall and compressibility effects, were uti- 
lized. The perturbational airloads were obtained utilizing perturbations of 
the three sectional aerodynamic coefficients with respect to both angle-of- 
attack and Mach number, as is outlined in a section of Reference 1. Compar- 
isons of Figure 107 with 109j and 108 with 110 show that the basic aeroelas- 
tic behavior identified by the eigensolutions with linear airloads is carried 
over and is but modified with the inclusion of perturbational airloads with 
stall. Generally, the effect of stalled airloads is to aggravate the stabil- 
ity bounds so identified; the same general trends with collective angle are 
qualitatively duplicated, but occur at lower collective angles. Again, these 
results are in good qualitative agreement with previous findings (Reference 
20), i.e., that- for either subcritical or supercritical configurations the 
flap or lag eigenvalue branch with the higher frequency is the one to go 
unstable. 

The following observations can be made of the above discussed nonvacuum 
eigensolutions : 

1. In light of (l) the consistency of the present findings with those 
obtained by other investigators for considerably simpler but basi- 
cally hingeless rotor systems, and (2) the fact that the eigen- 
solution results generated herein by the analysis (Reference 1) do 
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not consider the higher order effects of the redundant analysis 
(which are inherent to the composite bearingless rotor, CBR, con- 
cept), it is most likely that the present findings are typical of 
all hingeless rotors and not limited to the CBR concept. 

2, The built-in twist distribution for this lotor together with the 
assumed variable inflow, tends to produce a relatively uniform 
angle of attack distribution, and since for a MCA 0012 airfoil 
the static stall angle tends to decrease with Mach number, the 
spanwise extent of stall should tend to propagate inward from the 
blade tip. Since Figures 109 and 110 show the instabilities 
occurring at collective angles above that wherein the 75 percent 
span station stalls, these instabilities represent conditions of 
deep stall which are not easily reached operationally. 

3. These results retain qualitative usefulness in the context of 
providing limited understanding of the aeroelastic behavior of 
rotor blades in forward flight . Momentarily high pitch angles 
and conditions near to and above stall can and do occur on the 

I retreating portions of the disk in forward flight. 

I 

Summary of Eigensolutlon Results - The principal findings gleaned from 
the above described eigensolutlon analyses are: 

1. The subcritical blade as originally configured is free of any n/rev 
resonance problems . 

2. The supercritical blade, when configured with a softened torque tube 
(25 percent nominal bending stiffnesses), is free of any n/rev 
resonance problems . 

3. The effects of torque tube bending stiffness and blade pitch angle are 
to produce highly complex coupled mode shapes with nonuniform phasing. 

4. The existence of hingeless rotor flap-lag Instability at high pitch 
angle and in deeply stalled hovering conditions, as reported in the 
literature, was confirmed. This hovering instability, in itself, 
is neither a result of the innovative features of the CBR nor poses 
any problem for the rotor within any realistic operational envelope. 
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Trimmed Forward Flight Cases 


The remainder of the aeroelastic analysis is confined to investigations 
of various of the rotor characteristics in forward flight. For these inves- 
tigations the eigensolution analysis is inadequate and the time-history solu- 
tion capability of the aeroelastic analysis program must be used. For such 
investigations, especially with the attendant required parameter variations, 
the rotor must first be simulated in (stable) trim configurations to have 
practicality. This subsection describes the results of the use of the time- 
history solution to obtain these trim configurations. In particular, the 
selection of the particular trim conditions is first discussed. Then the 
resulting blade bending moments, stresses and torsion moments experienced by 
the two basic blade configurations in trimmed straight and level forward 
flight are discussed in turn. 

Selection of Trimmed Flight Conditions - For present purposes a high 
speed condition (m. = 0.44) with moderate blade loading = 0.057) was 

selected. More specifically, the rotor was trimmed to a lift and propulsive 
force of 59610 N (13,400 lb) and 7430 N (167O lb), respectively, with hub 
moments less than 5420 Nm (4000 lb -ft) in magnitude; a high forward flight 
speed of 86.9 m/sec (169 kts) at standard sea level conditions was selected. 
Uniform inflow was assumed and the aeroelastic analysis was used in the auto- 
matic trim mode to adjust the three control angles and the inflow ratio iter- 
atively to achieve the above defined flight condition. Three basic trim con- 
ditions were calculated: one each for the subcritical and supercritical con- 

figurations discussed in an above subsection (case Nos. 1 and 2, respectively) 
and an additional case (No. 3) for the subcritical configuration. For this 
latter subcritical configuration, the push-rod attachment point was reposi- 
tioned further outboard to investigate the effects of increased pitch-flat 
coupling. These three basic flight conditions are summarized in Table X. 

Blade Moment and Stress Distributions - Figures 111 through 113 present 
the calculated distributions of blade one-half peak-to-peak ( H FTP) and 
median bending and torsion moments over the (total) blade span for each of 
the two basic blade configurations. In these figures the moment values pre- 
sented for the structurally redundant, flexbeam-torque tube spanwise stations 
(those inboard of the juncture station, f ^ 0.29) are those appropriate only 
to the respective primary elastic structures. Thus, in Figures 111 and 112 
the bending moments presented for stations inboard of the juncture are those 
experienced by the flexbeam; conversely, in Figure 113 the torsion moments 
presented for these inboard stations are those experienced by the torque tube. 
Explicit calculations of the complementary bending and torsion moment distri- 
butions (in the torque tube and flexbeam, respectively) were not available. 
However, their presence and the relative significance of structural redundancy 
is indicated by the jump magnitudes in the various moment distributions at the 
juncture station. The effect of structural redundancy for these configurations 
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is thus seen to be small for flatwise bending, but moderate for edgewise 
bending and for torsion. Not unexpectedly, the maximum blade bending moments 
are seen to occur at the blade root. The negative median flatwise bending 
moments near the blade tip shown in Figure 111 are due mostly to the inclu- 
sion of an 11.3 kilogram tip weight and a greater than optimal built-in pre- 
coning angle. 

Significant differences are shown in the respective bending and torsion 
moments for the subcritlcal and supercritical configurations. Figures 111 and 
1 1-2 show the supercritical configuration to exhibit both flatwise and edgewise 
bending moments which are generally lower than those for the subcritical conr- 
figurations. Of particular significance are the marked reductions in edgewise 
bending moments for the supercritical configuration, both with -g-PTF values over 
the entire span and with median values over the outer portions of the blade. 
These bending moment reductions are well recognized to be inherent in super- 
critical blade designs and result from the favorable "vibration absorber" and 
equivalent hinge characteristics of this type of design. Figure II 3 compares 
the torsion moment distributions obtained for the two designs. It shows that 
the torsional moment carried by the flexbeam for the supercritical configiora- 
tion is less than that for the subcritlcal as indicated by the reduced jvimp 
values at the juncture. This result is consistent with the reduction in non- 
linear ^EI torsion coupling obtained when approaching a matched stiffness 
design. Figure 113 also shows a generally outward shift of the spanwise peak 
of the torsion moment (disregarding the j’umps at the j’uncture) with the super- 
critical design; the torsion moment distribution for this design is thus more 
uniform over the span. 

The flatwise and edgewise stress distributions over the flexbeam span 
for the two blade configurations are shown in Figures ll4 and 115 . The curves 
shown in these figures are obtained from the bending moment distributions of 
Figures m and 112 divided by appropriate l/c values. The most remarkable 
result shown in these figures is that in contrast to the moment distributions 
the stress distributions for the two blade configurations are quite similar. 
Indeed, the vibratory (^PTP) edgewise stresses for the supercritical design 
are even somewhat larger than those for the subcritlcal, while the median 
edgewise stresses for the supercritical are significantly larger than the 
median stresses for the subcritical. From the results of the above Composite 
Materials Investigation and Pull Scale Design sections, maximum allowable 
bending and torsion stresses of 379 Mw/m^ (55jOOO psi) and 24 MW/m^ (3500 psi), 
respectively, were established. Figures U4 and II 5 thus show that the blade 
designs are satisfactory from a root bending stress standpoint. 

Azimuthal Time Histories of Critical Stresses and Loads - From Figures 
ll4 and 115 it can be seen that the critical bending stress point is near the 
blade root and is herein taken to be the center of the first blade segment: r = 
0 . 054 . Pigiire 116 compares the azimuthal variations of the flatwise and edgewise 


50 


stresses at the blade root for both the subcritical and supercritical 
configurations. It can be noted that the flatwise stresses for the two con- 
figurations are quite similar and, commensurate with zero trimmed hub moments, 
show predominately 2P and 3P content. Comparison of the edgewise stresses for 
the two configurations shows predominantly IP content for each configuration 
but with opposite phasing. Since the edgewise bending dynamic system is lightly 
damped, a change of phase angle of nearly l80° is to be ejcpected for responses 
whose natural frequencies are either greater than or less than predominant IP 
forcing frequency. 

An additional stress item of interest is the critical flexbeam torsion 
(shear) stress; this stress is t 3 q>ically located at the blade station where 
the twist rate is maximum which occurs where the sectional dimensions are 
minimum. Since the flexbeam is tapered over its entire length, this station 
is herein taken to be the outer boundary of the outermost segment of the flex- 
beam; f s 0.288. Althotigh explicit torsional stress distributions over the 
flexbeam span were not available from the analysis, the torsion moment in the 
flexbeam immediately inboard of the junctTore was available as part of the cal- 
culation of the redundant analysis. This torsion moment was converted to a 
maximum torsion stress within the section using the sectional dimensions and 
standard stress analysis formulae (Reference 21). Figure 117 compares the 
flexbeam torsion stresses at the jiincture for the two blade configurations. 

These stresses are seen to be almost entirely IP and reflect the predominantly 
Bi cyclic control input. The relative lack of bias in these flexbeam torsion 
stresses is due to the assumed flexbeam pretwist of 13.5 deg (corresponding to 
a b^^ilt-in collective angle of 11.66 deg). The 
flexbeam torsion stress for the subcritical and 

and 22.6 Mw/m^, respectively. Although below the allowable torsion stress of 
24 MS/isl, these values do not leave the margin initially predicted in the Pull 
Scale Design section. The reason for this discrepancy is that the preliminary 
design stress calcTilation was based on the assimaption of a constant section 
flexbeam. The actual design, however, does have substantial taper giving rise 
to a significant increase in torsion stress at the outboard end: In this case 

that increase is approximately 23 percent. 

In addition to the critical stresses described above, an important dynamic 
response variable is the p\ish-rod load. Besides reacting those loads required 
to twist the flexbeam cyclically, the push-rod must equilibrate the blade in 
torsion and, hence, react the inertial and aerodynamic pitching moment distri- 
butions on the blade. The azimuthal variations in push-rod load, shown in 
Pigvire 118 for the two blade configurations, reflect these load characteristics. 
Comparison of Figure U8 with Figure 117 shows the steady and IP content of 
the push-rod load to be closely phased with the steady and IP torsion moment 
required to twist the flexbeam to achieve the trimmed control angles. Note 
that, due to the trailing edge push-rod attachment location selected for these 


half peak-to-peak values of 
supercritical are 22.8 MN/m^ 
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designs, a positive (up) push-rod load produces a negative (leading edge down) 
flexheam torsion stress. In addition, however, significant higher harmonic 
content is evident in these push-rod load time-histories which result from the 
outer blade inertial and aerodynamic moment loadings. Of particular interest 
is the intermittent (approximately) lOP excitation of the torsion mode in the 
retreating blade quadrants. 

Effects of Simplified Modeling of the Structinral Redundancy - As was stated 
above, the nominal calculation of the various time-history solutions comprising 
the forward flight portions or the aeroelastic analysis assumes the complete 
modeling of the structural redundancy (see Reference 1). Briefly the struc- 
tiiral redundancy arises from the presence of dual load paths, over the flex- 
beam-torque tube span, for both bending and torsion loads. Ihe redundant 
analysis formulated for the solution of this structural problem involves a 
full interaction of all bending and torsion loads at the juncture point as 
well as a statement of consistent torsion deformation at that point. An 
alternate simplified formulation of the structural dynamics of this portion of 
the rotor blade is afforded by making the following approximations : 

1. ■ The torque tube consists of an infinitely rigid massless link, hinge 

connected to the span at the juncture by equivalent rotary springs in 
the flatwise and edgewise senses. 

2. The flexbeam is torsionally equivalent to a simple linear rotary 
spring. 

3. The nonlinear torsion excitation (due to AEI bending-deflection 
effects) over the flexbeam span is approximated by weighting the 
usual nonlinear elastic torsion load distribution by the pseudo- 
torsion mode shape (see Reference 1). Since the pseudo-torsion 
mode is a monotonically Increasing function less than vinity over 
the flexbeam span, this weighting results in an attenuation of this 
term. 

It is to be noted that the Incorporation of this alternate simplified formula- 
tion reduces the analysis in large measure to. that for a "conventional” hinge- 
less rotor. Remaining differences are the time variation in structural twist 
and the attenuation of the nonlinear torsion excitation for the composite 
bearingless rotor. 

Using these nonredundant analysis simplifications together with the 
control angles defined in Table X for trim cases 1 and 2, two additional time- 
history cases were generated. These cases are herein denoted lA and 2A for 
the subcritical and supercritical configurations, respectively. It was found 
that substantially the same performance was achieved for cases lA and 2A as 
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for cases 1 and 2; the nonredundant analysis cases were equally found to 
satisfy the trim criteria enunciated above. The azimuthal variations of the 
critical stresses for the nonredundant analysis are compared with those ceilcu- 
lated using the complete redundant analysis in Figure 119a. This figure shows 
negligible differences in the flatwise and edgewise bending stresses and in 
the flexbeam torsion stresses for the two calculation methods. Substantial 
differences are noted, however, in the azimuthal variations of the push-rod 
load shown in Figure 119ti. It is significant to note that the respective 
azimuthal variations, with and without the full redundant analysis, show simi- 
lar higher harmonic signatures and differ mainly in their lower harmonic 
contents. 


Aerodynamic Stability Derivatives 

Two basic requirements exist for obtaining partial derivatives of the hub 
loads with respect to the control angles and components of hub velocity: First, 

such derivatives are needed for determining the static stability characteris- 
tics of the entire aircraft, as well as the control power characteristics of 
the rotor. Second, the derivatives can be used to generate trim derivative 
matrices for obtaining the required control angles and inflow for perturbations 
in the flight conditions. The following subsections discuss the aerodynamic 
stability derivative matrices defined by these requirements; results are 
obtained for the three basic trim cases discussed in the previous section and 
summarized in Table X. 

Basic Static Stability Derivative Matrices - For each of the three basic 
trim cases a matrix of partial derivatives is obtained; the partial derivatives 
are numerically obtained by perturbing (both positive and negative perturba- 
tions) the control parameter inputs of each trim case using the time-history 
solution of the analysis. The partial derivatives of each of the six (6) hub 
loads (3 forces and 3 moments) with respect to the three (3) control angles, 

, ^ , and 0 inflow ratio and advance ratio form a 6 x 5 matrix of 

fundamental derivatives. These derivatives are then modified to reflect the 
substitution of shaft angle and (nondimensional) forward flight speed as the 
fourth and fifth Independent control parameters. The results of perfoiming 
the numerical differentiation and modification are presented in Table XI 
where the various static stability derivative matrices are defined as: 



Pj a (A]^^, 0.75|J> ffg, V) (angles in degrees) 
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Comparison of the static stability derivative matrices of trim cases 1 
and 3 shows that the net effect of increasing the push-rod attachment point 
radius (to achieve increased pitch-flat coupling for the first flatwise mode) 
was actually to reduce the control power of the rotor. This can be seen by 
comparing elements Si2> 3^2 j S 41 for these two cases. This con- 

trol power reduction is contrary to what one would typically expect consider- 
ing only the increased (positive) pitch- flat coupling of the first flatwise 
mode. This behavior can be explained. by noting in Table X for case 3> the 
significantly increased (negative) pitch- flat coupling of the second flatwise 
mode associated with the outwardly repositioned push-rod attachment point. 
Since for trim conditions 1 and 3 the first harmonic responses of the first 
and second flatwise modal responses were, in fact, calculated to be in phase, 
the pitch-flat coupling of the second flatwise mode is thus seen to be pre- 
empting that of the first mode. 


Trim Derivative Matrices - For purposes of calculating additional trim 
cases a portion of the above defined fundamental derivative matrices for the 
three trim cases were inverted to form the trim derivative matrices presented 
in Table XII. These trim derivative matrices are defined as: 
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Transient Responses Due to Control Inputs 

Trim cases 1 and 2, whose parameters are given in Table X and whose 
stress characteristics are discussed in an above subsection, form the bases 
for investigating the transient blade responses due to control angle inputs. 

The effects of the control angle inputs were studied by combining the respec- 
tive trim control angles with separate time variable increments in each of the 
three control angles, 0 Aj_ and Bj_g. Those initial conditions on the 
response variables, appropriate to the unpertixrbed trim cases, were selected 
so that the deviations of the transient (stress) responses from the trimmed 
responses would provide direct measures of the effects of the incremental con- 
trol inputs. For each control angle two increments were made of the types 
shown in' Fig. 120. Both the ramp-step and pulse incremental inputs shown in 
this figure have a maximum amplitude of 0.5 deg with finite rise times to 
that amplitude. The responses due to the ramp-step input would be expected 
to be more severe as they represent off-trim conditions and would not generally 
converge to the original azimuthal stress response signatxires. 



The transient stress and co;itrol load responses due to these incremental 
J 3 ontrol angles are presented in Figures 121 thru 126 . For each of these 
figures, duplicate reference (trim) response time- histories are given for two 
periods (rotor revolutions comprising 36 O deg of azimuth). The initial tran- 
sient stress responses are given for these two revolutions and then the con- 
tinuing responses at subsequent revolutions are overlayed on these two revolu- 
tions, as indicated. Within any one plot the transient responses are plotted 
only up to and including the maximum (and mini mum) amplitudes. For those 
transients wherein several revolutions are required for convergence, only 
selected portions of the time-histories which include local, maxima are included. 
Within each of the following subsections describing the effects of each of 
the incremental control inputs, restilts are presented, in turn, for the sub- 
critical and supercritical configurations (trim cases 1 and 2 , respectively). 

Incremental Collective Angle, ^9.7^R - The results of the incremental 
ramp- step and pulse collective angle inputs are shown in Figxires 121a through 
122b. Observations to be made of these results are as follows; 

1. The responses all tend to show higher harmonic content than their 
respective trim responses during the initial transient revolutions. Thus, the 
incremental half peak-to-peak responses will not, in general, be linearly 
dependent functions of incremental input amplitude. 

2. With regard to the flatwise and edgewise bending stresses for both 
rotor configurations (Figures 121a and 122a), the converged maximum and mini- 
mum peaks are not attained until some elasped time (3 or 4 rotor revolutions) 
after the incremental collective input has been completed. The convergence 
tends to occur somewhat sooner with the subcrltical configuration, however. 

3 . Comparison of the flatwise bending responses for the two rotor 
configurations shows the transient maximum and minimum peaks for the subcriti- 
cal configuration to be significantly more negative than those for the super- 
critical configuration although the peak-to-peak amplitudes are comparable. 

4. Congjarison of the edgewise bending responses shows the supercritical 
rotor to be significantly more sensitive in edgewise bending to the incremental 
collective control angle input than is the subcrltical. 

5 . As shown in Figures 121b and 122b, both rotor configurations show 
negligible incremental flexbeam torsion stresses due to incremental collective 
angles other than that due to the additional 0.5 deg of flexbeam twisting. 
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6. Comparison of the transient push-rod loads for the two rotor 
configtirations shows the subcritical rotor to he significantly more sensitive 
in torque tube torsion (of which push-rod load is a measure) to the incremental 
collective control angle than is the supercritical. In particular, the 
increased sensitivity is confined to the retreating blade quadrants wherein 
the local incipient (9P - lOP) torsion excitation, already present in the 
trimmed responses, is substantially amplified. While a complete investigation 
of this phenomenon is beyond the scope of this study, it appears, from an 
examination of the available moment load distributions, that this increased 
higher harmonic excitation arises from the AEI nonlinear bending- deflection 
torsion excitation. As would be expected, the supercritical configuration 
which is more nearly "stiffness-matched" is less prone to this ts^e of excita- 
tion and it concurrently shows less sensitivity. 

Incremental Longitudinal Cyclic Angle, AAt ^ - The similar results of 
incremental ramp-step and pulse longitudinal cyclic angle inputs are shown in 
Figures 123a through 124b. As with the transient responses described above, 
those due to incremental cyclic angle also show higher harmonic content 
during the initial transient revolutions, than do their respective trim 
responses. Contrary to the incremental collective angle responses, those due 
to incremental A^^ cyclic angle are relatively benign and show little varia- 
tion from the trimmed values. Two main reasons exist for this subdued behavior. 
First, ■ incremental longitudinal cyclic for rotors with low flatwise stiffness 
(first flatwise modal frequency < l.lP) will produce mainly incremental first 
sine harmonic flapping. This flapping does not increase the effective ram 
inflow as does first cosine harmonic flapping which would result from both 
collective and cyclic angle increments. Secondly, because the Aj (and 
) control angles are impressed cyclically, the average incremental blade 
pifch angle around the azimuth due to incremental %s would inherently be no 
more than half that due to incremental collective. Consistent with previous 
results the effect of incremental A]_ control angle on transient push-rod 
loads are greater for the subcritical than for the supercritical. 

Bicremental L ateral Cyclic Angle_, ABj^, _ Figures 125a through 126b present 
the results of incremental ramp-step and pulse lateral cyclic angle inputs. 
Comparison of Figures 125a and 126a shows that the transient flatwise stresses 
due to the Bj_ inputs for the two rotor configurations are quite similar and 
are generally well converged within three rotor revolutions. Comparison of 
the transient edgewise stresses shows again that the higher harmonic content 
is increased with the application of the incremental control input; this 
characteristic is seen to be more prevalent with the subcritical rotor con- 
figuration. Comparison of Figures 125b and 126b shows that the flexbeam tor- 
sion stresses are little affected by the incremental B^ cyclic angle. Con- 
sistent with previous results, the transient push-rod loads for the subcritical 
rotor configuration are seen to be significantly more affected by the 
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incremental control angle than is the supercritical configiiration. It is 
to he noted that, for the subcritical configuration, the increased higher 
harmonic edgewise stress responses occxir over the same portion of the azi- 
muth (retreating blade quadrants) as do the push-rod load responses. 



Transient Characteristics Due to Partial Material Failure 


As is discussed in the above Composite Material Investigation section, 
composite materials typically fail by developing noncatas trophic cracks which 
reduce their elastic moduli. Thus, flexbeam failure becomes more of a struc- 
tural dynamics problem. That is, assuming a material failure (which will gen- 
erally be of the form of partial reductions in stiffness), can the structvire 
withstand the dynamic loads without destructive stress magnifications and/or 
aeroelastic instabilities? An assessment of the material failure problem was 
made by investigating the stress responses of the two blade configurations due 
to instantaneous reductions in the flexbeam bending stiffnesses. 

More specifically, the two blade configurations were first analyzed for 
their frequencies and mode shapes with separate 15 percent reductions in flat- 
wise and edgewise bending stiffnesses. The results of these reductions on the 
blade uncoupled natural frequencies are shown in Table XIII. Using these 
uncoupled nattiral frequencies and their associated mode shapes, together with 
the control angles and modal initial conditions for trim cases 1 and 2, tran- 
sient time-history solutions were obtained in a manner similar to the transients 
due to control inputs. In this case, however, the sources of input disturbance 
were the modal initial conditions. For lack of a more comprehensive stress 
analysis of failed composite materials the "failed" stresses were computed 
using the force-integrated bending moments together with the "unfailed" l/c 
values . 

Figures 127 and 128 present the transient flatwise and edgewise bending 
stresses for the subcritical and supercritical configurations, respectively. 

In each figure are shown the reference (trim) stress time-histories, previously 
discussed, together with the transient stresses accruing from instantaneous 
partial reductions in flexbeam flatwise or edgewise bending stiffness. As 
these figures show, the effects of these instantaneous stiffness reductions 
are practically negligible. Indeed, the flexbeam torsion stress and push-rod 
load time-histories were not included since the off nominal results could not 
be distingulshably plotted together with the reference time-histories. That 
so little incremental dynamic stress is generated by the instantaneous partial 
failures is consistent with the results of Table XIII. Therein it is seen that 
the stiffness reduction caused only small excursions in modal natural frequen- 
cies which were still reasonably well separated from integral harmonics. The 
small deviations of the stresses from the reference (nominal stiffness) 
responses can generally be explained on the basis of the changes in resultant 
dynamic magnification factors caused by frequency change. Most notable are 
the changes in edgewise stress for the two configizrations upon reduction of 
the flexbeam edgewise stiffness. Since the edgewise stiffness reductions 
brought the first edgewise natural frequencies closer to and further from IP, 
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respectively, for the subcritical arid supercritical configurations, it is 
reasonable that the IP edgewise stresses sho\d.d increase and decrease, respec- 
tively, for these two configurations. 


Effects of Uristeady Airloads 

An important aspect of the aeroelastic analysis of any helicopter rotor 
blade is the effect that xmsteady airloads have on the structural responses. 

Hie formulation of unsteady airloads used in the present aeroelastic analysis 
Is the variant of the (a. A, B) method described in Reference 22 . Briefly, 
this formulation states that the unsteady aerodynamic section coefficients are, 
for any fixed Mach number, generally nonlinear functions not only of angle-of- 
attack, a, but of nondlmensiontJ. angle-of-attack rate, A(=» c^/ 2 U), and angle- 
of -attack acceleration, B(» c 2 a/ 4 u^). While it has been noted in Reference 
23 that a proper simulation of stall-flutter requires the inclusion of vari- 
able inflow, it was decided that the effect of unsteady airloads should be 
determined using only constant inflow for a more meaningful compairison with 
the herein nominal quasi- static results. 

The application of unsteady airloads in the aeroelastic analysis was 
confined to an examination of its effects, with variations in blade loading, 
on the response characteristics of the subcritical configuration. Hie varia- 
tions in blade loading were achieved by straight-forward variations in the 
(constant) inflow, x, or equivalently, the rotor shaft angle-of-attack, ag. A 
meaningful comparison of the results obtained with and without unsteady airfoil 
is complicated by the fact that the various pertinent structural response 
quantities (^PTP stresses, etc.) are strongly dependent upon the state of rotor 
trim. Ihe approach selected for constnicting a meaningful comparison was, 
therefore, to make inflow variations on cases which were initially trimmed to 
the same hub loads. Accordingly, an additional trim case was calculated 
wherein the analysis included the unsteady airloads formulation; the control 
parameters found to trim the an8j.ysis with this incl^^sion are as follows; 

=» -1.22 deg, B^g » 9.39 deg, 9 =■ 12.09 deg, \ » -.0733). 

In Figure 129 are presented the results of the variations of rotor 
performance with shaft angle, oTg, using the aeroelastic analysis with and with- 
out the unsteady airloads formulation. In this figure are shown the two trim 
points about which the inflow variations were made. Note that, for comparison, 
an isolated point on each curve is presented wherein the control parameters 
appropriate to the norminal quasi- static airloads (case 1 ) were used in the 
analysis with the incltision of the unsteady airloads. A comparison of the 
Ci/ct curves with and without the unsteady airloads shows that, at the higher 
blade loading conditions, the performance calculated with unsteady airloads is 
rising at a greater rate and would appear to be “stalling out" at a higher 
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overall level. Ihe rotor drag levels with unsteady airloads are, furthermore, 
seen to be less than those with the quasi-static formulation. These res\ilts 
reflect the generally acknowledged effect of unsteady airloads on rotors: 
that dynamically stalled airloads create stall delays and hence lift over- 
shoots. These stall delays cause the airfoil sections to generally operate 
more efficiently than quasi-static theory wouild indicate in creating more 
average lift per (effective) rotor angle-of-attack. Further observations from 
Figure 129 are that the two trim points occur at different shaft angles and 
that the quasi- static trim point lies on a portion of the QjJa - o!q ciurve 
which has discernable curvatiire. From these observations it would appear 
that, as calculated, basic trim case No. 1 is a partially stalled rotor oper- 
ation condition. It should be noted that the higher power absorbed by the 
rotor with unsteady airloads is commensurate with its ability to pull more 
lift. 


Comparisons of the concurrent vibrational stresses and push-rod loads for 
the analysis with and without the unsteady airloads are shown in Figure 130. 

Of particular interest is the substantially higher vibratory flatwise stresses 
for the quasi- static formulation at the higher blade loadings. AlthoTigh much 
of this difference is attributable to the quasi- static results being further 
out of .moment trim than those for the unsteady airloads, it was nevertheless 
found that all harmonics of this stress through the tenth were substantially 
greater for the quasi- static than for the unsteady airloads formulation 
results. The comparison of ^PTP push rod loads in this figure shows a similar 
result for the two formulations. Here too, the results using the quasi-static 
formulation have substantially higher harmonic content through the tenth har- 
monic than do those using the unsteady airloads formulation. 

It is significant to note that these vibratory stress and push-rod load 
results are contrary to what one would expect based on established stall 
flutter theory. Stall-flutter is generally attributed to the energy extrac- 
tion mechanism provided by the hysteretic moment coefficient characteristic 
inherent in the unsteady airloads formulation. As such, stall flutter is 
generally a single degree-of- freedom torsional response phenomenon. The 
present findings indicate, however, that the onset of the excessive torsion 
responses and amplified flatwise bending responses found using quasi-static 
airloads is delayed with the incorporation of unsteady airloads. As was shown 
in an above subsection describing the effects of structural redundancy, at 
high frequencies, the CBR tends to behave much like a "conventional" hingeless 
rotor. It would thus appear that the excessive torsion and flatwise bending 
responses found at the high blade loadings, and previously found to a limited 
degree for incremental blade collective angle, constitute a potentially new 
phenomenon associated with the analysis of hingeless rotors using quasi-static 
airloads. As was shown in Reference 23 > a satisfactory calculation of stall 
flutter generally requires the inclusion of variable inflow which was lacking 
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in the present calculations. Thus, that stall-flutter per se was not 
encountered using the \msteady airloads is not prima facie evidence that stall- 
flutter would not occur on the composite bearingless rotor at these elevated 
blade loadings. An adequate exploration of this phenomenon is, however, beyond 
the scope of this study. 


Forward Flight Stability Analysis 

The aeroelastic stability of rotor blades operating in forward flight 
conditions is, as best, difficult to analyze for the basic reason that the 
governing differential equations of motion are significantly more complicated 
than those for which the usual mathematical tools for stability analysis were 
developed. In contrast to the familiar linear, homogeneous equations of 
motion with constant coefficients, those for the rotor blade in forward flight 
are highly nonlinear, contain periodic coefficients in various dominant linear, 
as well as nonlinear, terms and are substantially excited in all lower harmon- 
ics. Consequently, for any time-history response of finite duration the tran- 
sient responses often defy categorization as to exponential growth and/or 
transient convergence to periodicity. 

Mathematically, the stability of such nonlinear systems can be defined 
in terms of the linear stability of small perturbations away from some nominal 
periodic solution. This solution is but a converged time-history of the orig- 
inal harmonically excited nonlinear equations. The differential equations for 
the perturbations, however, are generally linear with periodic coefficients, 
and as such, constitute a problem in Floquet analysis for which several solu- 
tion methods are available. However, two major difficulties in applying this 
general technique to the differential equations of the rotor blade are as 
follows ; 

1. A converged time-history solution of the original nonlinear equations 
may not be available. Indeed, for cases of most interest, wherein 
the degree of stability is marginal, the (numerical) time-history 
solutions will not generally converge to periodicity within a prac- 
tical period of time. 

2. The linear perturbational eqmtions, to which Floquet theory must be 
applied, are not generally available in explicit form owing to the 
complexity of the original differential equations and to the afore- 
mention lack of converged solution. 

Becaoise of these two dlffictilties the stability of the rotor was analyzed 
herein using a combination of three basic ideas and/or mathematical tools; 
First, and most elementary, is to observe the time histories for responses 
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which clearly increase mondtonically in amplitude until unreasonably large 
amplitudes are reached. Time-histories of such responses could clearly be 
defined as unstable in the global sense, but quantification of the level of 
instability often becomes meaningless or undefinable. The second mathematical 
tool is to use the "frozen azimuth" eigensolution available in the aeroelastic 
analysis program to ascertain general coupled frequencies and trends in damping 
due to parameter variation. This second tool is used primarily as a gtiide for 
interpreting the results of the third mathematical tool available. This last 
and potentially most powerful mathematical tool for analyzing the rotor stabil- 
ity is the transient spectral stability analysis (TSSA) described in detail in 
Reference 24 and briefly in Reference 1. The results obtained from the TSSA 
are the values of the predominant frequencies inherent in the given time- 
histories as well as the eqiiivalent characteristic esjponents associated with 
each of these identified predominant frequencies. These equivalent character- 
istic exponents are interpreted in the same way as are the real parts of the. 
familiar eigenvalues defined for linear systems with constant coefficients. 

Two separate bodies of stability analysis results were calculated and are 
presented in this section. The first is an exploratory study of possible blade 
instabilities of one of the test models whose characteristics are described in 
previous sections. The second is an investigation of the transient stability 
characteristics of the full-scale design to complete the aeroelastic analysis 
section of this study. The model results are included in this subsection 
rather than in the model correlations section because they are entirely of an 
analsrfclc nature with a lack of quantitative data with which to correlate. 

These results furthermore serve as an introduction to and complement the simi- 
lar results obtained for the full-scale design. 

Exploratory Study of Model Rotor with Cantilevered Torque Tube - Of the two 
basic model rotor configurations tested, that with the cantilevered torque 
tube was selected as the one most likely to show adverse stability character- 
istics. This selection was made based on this configuration’s inherently 
softer torsion system arising from flexbeam bending effects (wobble mode), and 
on its pronounced pitch- flat coupling effects. Furthermore, it was found 
experimentally that the cantilevered torque tube configuration manifested more 
unusual aeroelastic response characteristics than did the pinned-pinned 
torque tube configuration, and would thus be more susceptible to potential 
instabilities. The remainder of this subsection describes stability results 
obtained for the cantilevered torque tube configuration. The pinned-pinned 
configuration was considered to be stable and, consequently wasn’t analyzed. 

In accordance with the above described approach for stability analysis, 
the aeroelastic analysis program was run in all three calculation modes: 
eigensolution, time-history solution, and transient spectral stability anal- 
ysis. Parametric variations were made on advance ratio and rotor speed for 
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two push-rod attachment positions- giving relatively large (+) and (-) values 
of -first mode pitch- flat coupling, respectively. Advance ratio, was varied 
up to a value of 0 . 4 ? at the nominal fixed value of rotor speed, fiR, of 99*1 
m/sec. Rotor speed -was vaxied up to the nominal -value holding advance ratio 
fixed at 0.4?; forward flight speed was thus varied in proportion to rotor 
speed to maintain fixed advance ratio. Leading and trailing edge push-rod 
attachment points were used which gave respective first mode pitch- flat 
couplings of 0.883 and - 0 . 654 , respectively. 

The eigensolution -was calculated at zero azimuth to more nearly obtain 
azimuthally average eigensolutions . All time-histories were calculated with 
zero control angles, an inflow ratio of -0.005, and (nominal) initial condi- 
tions of .02 on the -first flatwise and edgewise bending mode deflections only 
with all other modal initial conditions zero. The time-histories were calcu- 
lated for time periods siifficiently long to obtain data string channels of 512 
points for use in the transient spectral stability analysis. The time-history 
quantities selected for these data string channels were out-of-plane and in- 
plane deflections (z^ and y^, respectively) and torsional deflection (9g), all 
as calcTolated at the blade tip. 

The eigenvalue results for selected nominal tip speed cases are shown in 
Table XIV. Presented are the real and imaginary parts of the eigenvalues 
together with indications of the predominant modal components comprising the 
coupled mode shapes. Note that the 2 T modal component in these results denotes 
the rigid body torsion mode res-ulting from flexbeam bending effects and is 
herein referred to as the "wobble mode". Of most significance in these results 
is the sensitivity of the eigenvalues to the initial conditions and to advance 
ratio. These res-ults are useful primarily in identifying modal response fre- 
quencies, u), and in revealing destabilizing trends. Most noteworthy of these 
results is the ability of blade deflections to destabilize one of the coupled 
second flat-wise bending modes and the ability of advance ratio to coalesce the 
first flat-wise and edge-wise bending modes. Since these eigensolutions are 
"frozen azimuth" res-ults they are not generally indicative of the effective 
stability over the whole rotor azimuth. For this, the time-history solution 
together with the transient spectral stability analysis must be employed. 

The effects of advance ratio variation on modal frequency and blade 
stability are shown in Figiires 131 and 132 for the trailing and leading edge 
push-rod configurations, respectively. Based on the eigensolution results 
the frequency range of Interest was deemed to be from 0.8p to 3 . 5 Pj resonant 
frequency identification and equivalent characteristic exponent calculation 
were therefore restricted to this frequency range. For both configTorations 
the results indicate little -variation in modal frequency and generally stable 
characteristics, as indicated by the uniformly negative values of equivalent 
characteristic exponent, o. In the reduction and interpretation of the 
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transient spectral stability results various of the identified frequencies 
were omitted from the presented results. Those frequencies which were either; 
integral harmonics with negligible damping (responses due to harmonic excita- 
tion) or of nonintegral harmonic frequencies but with the same noninteger part 
and same characteristic exponent as other identified frequencies (image or 
nonfundamental resonances), were so omitted in the interests of clarity. 

Detailed comparison of Figures I 31 and 132 shows first that although both 
configurations are stable, the leading-edge attachment configuration (with 
positive pitch-flat coupling) has significantly less margin than does the 
trailing edge configuration. Also shown are the generally lower modal fre- 
quencies identified for the leading-edge attachment configuration. This result 
is consistent with known effects of positive pitch- flat coupling. Finally, it 
should be noted that although a single frequency of approximately 1.36 p was 
identified for each of the three channels, consistently different damping 
levels are indicated for the flapping (z^) channels as shown with dashed lines 
in the figures. It can fiurthermore be noted that these I. 36 P characteristic 
exponents (for the flapping channel) are more and less stable, respectively, 
for the trailing and leading edge push rod attachment configurations, than 
their corresponding inplane and torsion I. 36 P characteristic exponents. 

The effects of rotor speed variation on the stability characteristics of 
the trailing edge push-rod configuration are shown in Figures 133 and 13^. 
Figure 133 shows the variations in the "frozen azimuth" eigensolutions with 
rotor speed variation. Of most significance in this figure is the coalescense 
of the wobble mode (2 t) with the second flatwise mode (2F) at 3.2P for the 
nominal rotor speed of 99*1 m/sec and the associated indicated instability. 

As demonstrated in Table XIV this unstable eigenvalue is deflection dependent 
and is not applicable to conditions around the entire rotor azimuth. Figure- 
134, on the other hand, presents the transient spectral stability resu2ts, but 
in the same format as Figure 133 for direct comparison. It is to be noted 
that several features of the frequency variations are well duplicated including 
the frequency coalescence near 3P, and the first flatwise and edgewise mode 
frequencies. Comparison of the damping levels shows remarkable similarities 
despite the effects of the nonlinearities and periodic coefficients. Signifi- 
cantly, the destabilizing trend of the low frequency Inplane (edgewise) mode 
is well duplicated in the two figures. In Figxxre 134, the dashed line denoting 
damping results for the low frequency flapping mode indicates a quantitatively 
undetermined connection, established qualitatively on the basis of the fre- 
quency resiilts. The lack of damping determinancy at 85.34 m/sec for this 
mode resulted from the erratic nonmonotonic behavior of the amplitude of this 
frequency content within the flapping channel time-history. These figures 
also show remarkable similarity in predicting a general loss of damping with 
the coalesence of the flatwise and torsion modes. It should be noted that 
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well damped modes leave little vibratory content to be identified by the 
transient spectral stability analysis. Hence the well-damped modes fotind by 
the eigensolution will generally be missing from the TSSA resiilts. 

Bie effects of rotor speed variation on the frequency and stability 
characteristics of the leading edge push-rod configuration are presented in 
Figure 135* The eigenvalue results for this case were found not to differ 
significantly from those for the trailing edge configuration shown in Figure 
133 said hence, were omitted herein. Drawing upon those results, however, the 
identified frequency plot was constructed. The dashed curve was constructed 
by subtracting a nondimensional frequency value of 2 from the second curve 
(oc^) recognizing that the («2 curve coiold possibly be an image resonance, and 
that the so constructed oag-S curve connects with the identified frequency at 
the nominal rotor speed to create a smooth curve exactly where one should 
exist. Examination of the characteristic exponent curves and comparison with 
Figure 134 shows a general relative lack of stability margin for the leading- 
edge push-rod attachment configuration. Indeed, at a tip speed of 30.5 m/sec, 
the time-history solution calculated cotild be unqul vocally temed unstable. 

Not only were the characteristic exponents for all identified modal frequencies 
found to be positive, but the responses achieved unreasonably large magnitudes 
and were terminated in the fourth rotor revolution. 

The principal findings of the model rotor exploratory study are first 
that the trailing-edge push-rod config\jration with negative (first mode) pitch- 
flat coupling is significantly more stable than the leading- edge push rod con- 
figuration. Second, whereas mstable responses were calculated for the latter 
configuration, none were calculated for the former. Third, no instabilities 
were calculated for either configuration at the nominal rotor speed at any 
advance ratio up to the maximum tested value of 0.4?. These findings are thus 
consistent with the generally stable experimental results fotond within the 
test envelope. Biese calculations serve the additional piup>ose of introducing 
the stability analysis techniques also used in the following subsection. 

Transient Stability of Full Scale Designs - In contrast to the above 
discussed cantilevered torque-tube configured model rotor, the full-scale 
design configurations with inboard snubber demonstrated generally more stable 
response characteristics. Once again, the evaluation of these stability 
characteristics req\iired the consideration of results produced by the three 
bEisic mathematical tools available. For this portion of the study parametric 
variations were made about each of the two basic configuration trim conditions. 
Those parameters so varied were equivalent blade mass ratio (air density), 
blade loeuiing (rotor shaft angle with fixed controls), advance ratio, push- 
rod attachment station (pitch-flat coupling), and torque tube flatwise bending 
stiffhess. For each parameter variation, the blade response initial condi- 
tions were selected to be the appropriate trim values but with increments of 
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0.02 added each to the first flatwise and first edgewise bending modal degrees- 
of- freedom. Additionally, for some of the increased air density cases, the 
trim parameters were proportionately reduced to preclude unrealistic rotor 
stall conditions from obscuring the transient stability characteristics at 
issue. 


Results of the stability analysis for the full-scale designs are shown in 
Table XV and in Figures I 36 thru I 38 . Table XV shows both rotor design trim 
conditions to be stable from a "frozen azimuth" eigensolution standpoint. 
Furthermore, all calculated time-histories required to form the above enumer- 
ated parametric variations were found to be stable; that is, in every case the 
response amplitudes were well bounded and devoid of any divergent character- 
istics. 

The results of using the transient spectral stability analysis (TSSA) on 
the calCTOlated time-histories are presented in Figures I 36 thru 138 . These 
figures show the various identified resonant frequencies together with the 
equivalent (linear) characteristic ejspionent associated with each such frequency. 
Figure I 36 shows the results for variations in equivalent blade mass ratio, 

^equ. in rotor shaft angle, a^. Figure 137 presents the effects of 

varjn.ng advance ratio at the nominal value of mass ratio of 10.45 as well as 
at the augmented value of 12.54 (1.2 x nominal). Lastly, the results of 
varying push-rod attachment station (to achieve indirect variation in pitch- 
flat coupling) are shown in Figure 138. 

Interpretation of the results of Figures 136 thru I 38 must be made not 
only in light of the generally stable indications afforded by the eigensolu- 
tions and by the stable boundedness of the time-histories, but also by con- 
sidering the selected operation of the transient spectral stability analysis. 
Typical TSSA operation Involved utilizing all of approximately l4^ azimuthal 
revolutions of time-history for frequency identification. However, the selec- 
tion of only 100 time shifts in the calculations for the characteristic expon- 
ents appears to have been too few because it effectively limited this calcula- 
tion to a relatively short window within the time histories from the start of 
the revolution 6 thru to 75 percent of revolution 8. Another factor to be 
considered in the interpretation of these results is that the TSSA identified 
significant Floquet-like behavior. Ey far the most consistently identified 
mode was the lightly damped first edgewise mode. However, as is shown in 
these figures, other modes were identified which had frequencies separated 
from the edgewise mode by integral harmonics and which were generally found to 
have approximately the same characteristic exponents. A general feature of 
these figures, however, is that most of these multiple or image roots demon- 
strate stable (negative) characteristic exponents. A significant point to be 
made is that although the basic supercritical trim case was calculated to 
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exhibit positive (unstable) characteristic exponents, it is, nonetheless, 
stable since it did in fact converge to periodicity. 

The principal findings of the full-scale design transient stability 
analysis are first that, for i parametric variations, both designs were 
found to be stable. This evaluation is based principally on the fact that all 
time-history responses were found to be bounded for time periods corresponding 
to sixteen (l6) rotor revolutions. Second, as a quantitative analysis tool, 
use of the transient spectral stability analysis to predict stability levels 
was generally unsuccessful. Due to an inproper selection of operational 
parameters the TSSA did not use sufficient data for resolution of the calcu- 
lated characteristic exponents. 

An addendum to these findings is that in the course of making these 
parametric variations it was found that some of the increased air density 
cases, which used the nominal trim parameters, were found to produce -unbounded 
responses. It is conjectured that these responses are stall related and simi- 
lar to others obtained under high blade loading conditions and discussed in 
the "Transient Responses due to Control Inputs" and "Effects of Unsteady Air- 
loads" subsections. The mechanism of this phenomenon is not fully imderstood 
and -was not explored herein since it -was beyond the scope of this study. 
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CONCLUSIONS AND RECOMMENDATIONS 


1. Based on model rotor wind tunnel experiments, the Composite Bearingless 
Rotor with pinned-pinned control torque tube, is stable over a tested 
range of flight conditions and exhibits performance and blade response 
characteristics similar to conventional hingeless rotors. 

2. Based on static and fatigue tests , graphite/epoxy is ideally suited to 
satisfy the strength and modulus requirements of the Composite 
Bearingless Rotor (CBR). 

3. A relatively low stiffness rigid body torsion mode is inherent in the 
cantilever torque tube design causing significant pitch washout and 
potentially destabilizing pitch-flap coupling. 

4. Based on a preliminary design study, a reduction in rotor system wei^t 
of 23 percent can be achieved with the CBR as compared to a conventional 
articulated rotor of the same size. 

5. Thfe unique structural characteristics of the CBR flexbeam, with its 
azimuthally varying twist distribution, can be analytically modelled. 

6. As predicted by theory, and confirmed by experiment the effect of flex- 
beam twist on flat to edge coupling for the pinned-pinned configurations is 
minimal. The edgewise responses appear to be dominated by the 1 P and 
2 P drag airloads. 

7. Use of uncoupled modes together with appropriate twist related displac- 
ment correction functions provides a viable method for obtaining 
coupled mode response characteristics of rotor blades with arbitrary non- 
linear (continuous or discontinuous) and time-variable twist as is 
typical of the composite bearingless rotor. 

8. The significance of the rigourous modeling of the structiiral redundancy 
features of the CBR concept appears to be limited to the lower harmonics 
of the torque tube torsion stresses and, consequently, the push-rod 
loads . 

9. Of the two torque tube configurated rotor models tested, the pinned- 
pinned and cantilevered, the former is the more practical for immediate 
further development. In contrast to the latter, it generally demonstrated 
typical hingless rotor characteristics consistently devoid of aeroelas- 
tic anomalies. 
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10. Although the cantilevered torque tube configurations of the model rotor 
were generally stable within the test envelope, the leading edge attach- 
ment configurations would appear to have less stability margin than the 
trailing edge configurations and could lose aeroelastic stability 
entirely at conditions of concorrrent reduced rotor speed and moderate 
advance ratio. 

11. Based on the results of the model rotor tests and preliminary attempts 

to achieve trimmed stable flight conditions for a full-scale cantilevered 
torque tube design, the alternate full-scale design incorporating a 
snubber torque tube can more readily and practically be developed. (from 
this basic conclusion it was further concluded that only the snubber 
design would be included within the scope of the full-scale study.) 

The following additional conclusions relate to the evaluation of the 
full-scale (snubber torque tube) designs on the basis of limited ^aero- 
elastic considerations: 

12. The use of the CBR concept in a full-scale helicopter main rotor appli- 
cation is practical. 

13. Neither of the two configurations considered, one subcritical and one 
supercritical, are susceptible to any instabilities by virtue of the 
unique CBR features. While various instabilities were demonstrated for 
extreme combinations of blade loadings and advance ratio, they appear to 
be phenomena generally characteristic of all hingeless rotor systems. 

14. Both full-scale designs can be configiured to be free of any significant 
n/rev coupled frequency resonances. The bending stiffening of the 
torque tube needs to be considered and carefully controlled for the 
supercritical (soft edgewise) design to avoid l/rev dynamic magnifica- 
tions. 

15. The selected trimmed flight conditions at high speed (p, = 0.44) and 
moderate blade loading {C-^/a = 0 . 057 ) represent partially stalled con- 
ditions and, hence, ones near to an operational stall boundary. 

16. The designs are practical from the standpoint of flexbeam vibratory 
bending and torsional shear stresses. The flexbeam torsional stresses 
are marginally satisfactory at the hi^ speed design condition, and an 
additional design iteration would be justified if greater margin were 
desired. 

17. Transient blade bending stresses induced by perturb at ional control 
angles and partial material failures (15^ reductions) are not critical . 


69 


At those combinations of blade loading and high forward speed defining 
the stall boundary blade torsion moments and, consequently, push-rod 
loads could become critical due to interactions of nonlinear AEI tor- 
sion excitation (accruing from hingeless rotor flatwise x edgewise bend- 
ing) together with unsteady aerodynamics (dynamic stall and variable 
inflow). Because of the more nearly "stiffness-matched" characteristics 
inherent in the supercritical design, this design should be less 
critical. 

Although manifesting lower vibratory edgewise bending moments , the 
supercritical design can exhibit vibratory edgewise stresses greater 
than or equal to those for the subcritical design. The supercritical 
design can, furthermore, demonstrate greater sensitivity to perturba- 
tional control angles. 

Nominal bending stiffnesses of the torque tube can be sufficiently high 
as to result in significant portions of the bending moments at the 
juncture being carried by the torque tube, and in significant modifica- 
tions of the coupled natural frequencies. 

The examination of time-history solutions for boundedness was the most 
reliable indicator of aeroelastic stability employed. The transient 
spectral stability analysis, which utilizes Fourier Analysis techniques 
together with the time-history solutions, produced inconsistent results. 
This lack of consistency is most probably due to an inefficient implemen- 
tation and subsequent utilization of this basic and potentially powerful 
concept. 

The following recommendations are suggested which would extend the 
research in bearingless rotor systems. 

Exploration into the feasibility of a "matched stiffness" bearingless 
rotor should be conducted. The rewards of further weight reductions and 
control loads may be achievable with such a design. Such a design is 
necessarily a highly supercritical design and the attendant potential for * 
air and ground resonance difficulties must be solved. The application of 
supplementary integral material damping may provide an important solution 
to this concept. 

The analysis should be extended to include (a) the calculation of torque 
tube bending stresses, and (b ) explicit loadings on the torque tube. In 
direct support of a "matched stiffness" design the analysis should 
furthermore be extended to include (c) the hub degrees -of-freedom with 



multi-blade capability, (d) a modelling of hysteretic structural damping 
to simulate the application 'of discrete elastromeric dampers or supple- 
mentary integral material damping, and (e) control system dynamics. 

Such an analysis would thus go beyond the capabilities of the eigensolu- 
tion and Floquet approaches to include the full influences of nonlinear- 
ities, unsteady stalled aerodynamics, and realistic material damping. 

3. Detailed studies of the destabilizing pitch-flap coupling mechanisms 
should be conducted to determine if an exceptable solution to the canti- 
levered torque tube design exists for main rotor applications. 

4. Indications of a potentially critical phenomenon of hingeless rotors have 
been observed in the calculations. This phenomenon is stall-flutter like 
in that it involves the higher frequency torsion mode and occurs at high 
speed near the stall boundary. In contrast, however, it was calculated 
to occur principally when using only quasi-static airloads, involves 
appreciable blade bending and, hence, significant torsion excitation 
through the nonlinear AEI term. Whether this phenomenon is an analytic 
anomaly or a real problem area of hingeless rotors should be resolved. 

5. Further development of the transient spectral stability analysis should 
be made to assess its full potential and to increase its reliability in 
routine usage. 
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TABLE I 


METRIC TO ENGLISH CONVERSION FACTORS 
Bv 


To Find 


39.370 



— .7376 

3L8.5 

— .22U8 

— .1U5O 


— Ib-ft 

— — — lb-in. 2 

— — — lb-sec2 

— — — (lb/in.2)xl0 


TABLE II 

Unidirectional HTS Carbon/ Epoxy Static Test Resiilts 


Material 

Lot 

Number 

lA-11 


lC-4 


lD -52 


lD -52 

lA 

lD -52 

IB 


Shear 
Strength 
k 



99.3 

(14.4) 

112.4 

(16.3) 

115.8 

(16.8) 

113.1 

(16.4) 


Shear 

Modulus 


6.48 

(0.94) 

6.55 

(0.95) 


Flexural 
Strength 
s 


1.5 

( 222 ) 

1.76 

( 255 ) 

1.53 

( 222 ) 

1.67 

(242) 


Flexural 

Moduliis. 


156 

( 22 . 6 ) 

151 

(21.9) 

152 

( 22 . 1 ) 

144 

(20.9) 

132 

( 19 . 1 ) 


* Not measured 
















TABLE III 


Combined Loading Fatigue Results 
HTS Graphite/Epoxy Specimens 


No. 

^edge 

C'flat 

axial 

^ total 


KC4-4 


KC4-3 


KC4-2 


KC4-1 


KCI-6 


KCl-5 


G-/Gi Ef 


stress 

units : 

MN/m^, 

(Ksi) 

• 




173.7 

(25.2) 

496 

(72.0) 

0 

670 

(97.2) 

30 

(4.3) 

107 

.96 

.85 

122.7 

(17.8 

436 

(63.2) 

0 

558. 

(81.0) 

26 

(3.8) 

107 

1.00 

1.00 

87.6 

(12.7) 

369 

(53.5 

0 

456 

(66.2) 

22 

(3.2) 

10^ 

.84 

1.00 

81.6 

(12.7) 

369 

(53.5) 

0 

456 

(66.2) 

22 

(3.2) 

107 

.94 

1.00 

173.7 

(25.2) 

496 

(72.0) 

69 

(10) 

739 

(107.2) 

30 

(4.3) 

16^ 

.51 

.98 

84.1 

(12.2) 

341 

(49.5) 

69 

(10) 

494 

(71.7) 

21 

(3.0) 

107 

1.00 

1.00 

135.8 

(19.7) 

382 

(55.4) 

i 

69 

(10) 

587 

(85.1) 

23 

(3.3) 

10'^ 

1.00 

1 

1.00 


964-4 


» Split at clamp due to excessive tightening in trying to prevent slipping 










TABLE IV 


MODEL SCALE BEAEINGLESS ROTOR DESIGN PARAMETERS 


Parameter 

Torque Tube Design* . . . ... • • • . 

Design Tip Speed, R, m/sec . . , . . 

Ijlo . Of Blades ,b 

Radius, B,m 

Chord, c, m*« 

Solidity, O' 

Blade Root Offset ... 

Preconing , a^ , deg 

Flexbeam Length . . . 

Airfoil* 

Linear Twist , 6 ^ deg 

c.g. Location* 

e . a . Locat ion* 

Tip Mach No. , M • 

Tip Reynolds No., Re 

Lock No. , Y 

Froude No., n^R/g 

1st Flapwise Natural Frequency, (5^^ . . 

2nd Flapwise Natural Frequency, . . 

3rd Flapwise Natural Frequency, . * 

1st Edgewise Natural Frequency, , . 

2nd Edgewise Natural Frequency, . , 

1st Torsional Natural Frequency, uJ . . 

^1 • 


Model Rotor Values 


Cantilever (Fig. 9) 

99.1 . . 

4,2 

0.610 . . 

0.0391 

O.OBlT, 0.0408 . . 

0.05R . . 

2 . . 

0283R . . 

NACA 0012 . . 

0 . . 

0.25c . . 

0.25c . . 

0.291 . . 

0.287 X 10^ 

4.27 . . 

1640. . . 

1.09 . . 

2.90 . . 

5.66 . . 

1.35 . . 

8.52 . . 

7.15 . . 


, Pinned (Fig. 8) 

. 99.1 
. 0.610 

.. 0.0391 

. 0.0817, o.o4o8 
. 0.05R 

. 2 

. 0.283R 
.. NACA 0012 

. 0 

. 0.25c 
. 0.25c 
. 0.291 

. 0,287 X 10^ 

. 4.15 

. 1640. 

. 1.09 
. 2.66 
. 4.70 
. 1.35 
. 7.72 

. 7.02 


* Outboard of Torque Tube 



TABLE V 


/ 


/ 


SCALE FACTOR SUMMARY 


Parameter 

Factor 

Weight or Mass per Unit of Length 

(Rfs/Rms)^ 


Total Mass 

(%sAms)^ 

R^ 

Moment of Inertia per Unit Length 

(Rj-s/Rms)^ 

R^ 

Total Inertia 

(%s/Rms)^ 

P 

Force 


r^5^ 

Moment 

(Rfs/Rms) (^R)fs^/(^^^ms^ 


Stress 

[(or) /(hr) 

PS MS 


Stiffness 

(Rps/R]yjg)^ (QR)^g2/(nR)^2 

-6-2 

R n 

Nondimensional Blade Natural 
Frequencies 
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TABLE VI 

FORWARD FLIGHT TEST COKDITIONS 
QR= 99-1 m/sec 


h 


Collective Angle, 
®.75 » 

Torque Tube Attachment 

2 

0.25, 0.35, 0.47 

4 

Pinned 

ti 


8 

It It 

k 

0.25, 0.35, 0.47 

4 

It tt 

It 

0.25 

8 

tt tt 

It 

0.35 

It 

tt tt 

ft 

0.25 

9.3 

tt It 

n 

0.35 

10.1 

It It 

tt 

0.47 

It 

It tt 

2 

0.25 

7.4 - 7.5 

Cantilever (t.e.)* @ r=0.060R 

tt 

0.35 

It IT 

It It 

tt 

0.47 

It tl 

It tt 

It 

0,25 

13.5 - 13.6 

It It 

It 

0.35 

It It 

tt It 

tt 

0.47 

Tt tt 

It tt 

It 

0.35 

16.6 - 17.5 

tt It 

4 

0.25 

7.5 

. ft II 

2 

0.25 

6.6 

Cantilever (t.e.) @ r=0.080R 

If 

0.25 

12.5 

It It 

4 

0.25 

7.4 

II ft 

2 

0.25 

6.1 

Cantilever (l.e.)* @ r=0.091R 

tt 

0.35 

6.5 

It ft 

ft 

0.47 

5.5 

It n 

tt 

0.25 

9.5 

It It 

tt 

0.25 

11.1 

It It 

tt 

’ 0.35 

11.1 

It IT 

It 

0.47 

11.1 

It tl 

tt 

0.35 

13.2 

It It 

tt 

0.47 

14.4 

It tt 

4 

0.25 

7.4 

It tt 

2 

0.25 

6.5 

Cantilever (l.e.) @ r=0.070R 

2 

0.25 

9.7 

It tt 


NOTE: At each test condition Ofg was generally varied over a range of eight or 

ten degrees in two degree increments, 

* t,e. - trailing edge; l.e. = leading edge 
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TABLE VII 


PARTIAL DERIVATIVES OP PERFORMANCE WZffl 
RESPECT TO CONTROL ANGLES 


(a) Lift Coefficient/Solidity, Cj/a 


Torque Tube Type 

M- 

6/60-3 


6/6A^3 

6/6B1S 



deg“^ 

Pinned-Pinned 

.25 

.00614 

.01710 

-.00111 

-.00674 


.35 

.00928 

.02161 

-.00167 

-.01038 

Cantilevered (t.e.) 

.25 

.00443 

.00934 

.00077 

-.00288 

@ r = , 060 R 

.35 

.00540 

.01029 

.00098 

-.00484 

Cantilvered (l.e.) 
@ r = .09IR 

.35 

.00874 

.01639 

-.00272 

-.00950 


(b) Drag Coefficient/Solidity, C^/a 


Pinned-Pinned 

.25 

-.00038 

-.00226 

.00040 

.00056 


.35 

-.00193 

-.00294 

.00056 

.00167 

Cantilevered (t.e.) 

.25 

.00029 

-.00091 

.00026 

.00004 

@ r = .06OR 

.35 

-.00021 

-.00119 

.00003 

.00057 

Cantilevered (l.e.) 
@ r = .09IR 

.35 

-.00130 

-.00225 

.00025 

.00122 


(c) Torque Coefficient/Solidity,. Cq/a 


Pinned-Pinned 

.25 

0 

-.00018 

.00048 

.00114 


.35 

.00036 

.00030 

.00048 

.00037 

Cantilevered (t.e.) 

.25 

.00029 

.00033 

.00010 

.00032 

@ r = .O6OR 

.35 

.00024 

.00026 

.00008 

.00023 

Cantilvered (l.e.) 

@ r = .09IR 

.35 

-.00048 

-.00066 

-.00031 

.00062 




TABLE VIII 


MODEL TEST PARAMETERS DEFINING CORRELATION CASES 


Group 

Case 

b 

Torque-Tube 

Configuration 

Geometric Pitch-Flat 
Coupling, 6 Vi , i =1, 2 ,3* 

Density 

Ratio 

Velocity 

km/sec 

Geom, Shaft 
angle, deg. 

®.75 

deg. 

^8 

deg. 

Bis 

deg^ 

I 

i 

f 25.03 

2 

pinned 

0., 0. , 0. 

.91*98 

88.96 

-6.0 

8.2 

-2.3 

5.3 


1 

87.08 

1 

+ 





.9571 

88.96 

-5.0 

7.9 

-2.1 

6.1 



22.03 







.9656 

88.96 

-8.0 

9.8 

-2.3 

5.1* 

II 


22.0lt 







.9638 

121*. 51* 

-8.0 

9.6 

-1.0 

6.1 



1 21.05 







.9603 

167.35 

0. 

3.1 

-1.9 

5.2 



22.05 







.9585 

167.72 

-8.0 

9.6 

-1.6 

6.2 

III 

( 11.05 







.9380 . 

124.90 

-6.0 

4.0 

-0.4 

1.4 


t 11.06 







.931*5 

125.09 

-6.0 

4.0 

-0.4 

2.7 



82.10 





cantilevered 

-.656,- 

^62^93*^^ 

.9562 

89. ll* 

0. 

6.7 

-0.4 

4.8 

IV 


83.06 






t 

.9579 

89. 14 

-8.0 

11.0 

0. 

6.5 



72.06 





+.422,1 

.302,-11.0*** 

.9253 

168.46 

0. 

5.6 

-0.2 

8.0 



73.08 






n 

.9187 

168.09 

-8.0 

14.5 

-0.7 

10.8 


^Geometric pitch-flat coupling is computed from flatwise modal distributions and torque tube geometry; cantilevered 
torque tube is attached to blade spar 19*05 cm from rotor center. 


**Push-rod/torque tube attachment point is 4.90 cm from rotor center and 2.97 cm behind feathering axis. 
***Push-rod/torque tube attachment point is 5.54 cm from rotor center and 2.96 cm ahead of feathering axis. 
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TABLE IX 


EXPERIMENTAL AND ANALYTIC ROTOR PERFORMANCE FOR CORRELATION CASES 


Case 





Cl/<^ 

Cjj/a 

Oo/» 

exp. 

anal.* 

exp. 

anal.* 

exp. 

anal.* 

25.03 

.248 

-.00374 

(-.0315) 

-.02518 

.00462 

.07832 

.07832 
( .04837) 

-.01079 

-.00877 

(-.00565) 

.00925 

.00643 

(.00532) 

87.08 

.249 

-.00373 

(-.0331) 

-.02232 

.02366 

.07410 

.07416 
( .00388) 

-.00898 

-.00735 

(-.00481) 

.00583 

.00606 
( .00501) 

22.03 

.247 

-.02754 

(-.0366) 

- .02658 

.00469 

.07079 

.07069 
( .05931) 

-.01137 

-.00992 

(-.00808) 

.00793 

.00748 

( .00675) 

22.0k 

.346 

-.03201 

(-.0483) 

.00087 

-.01807 

.05027 

.04931 

(.04166) 

-.00609 

-.00549 

(-.00419) 

.00573 

.00693 

(.00608) 

21.05 

.469 

.03275 

(.0064) 

-.02248 

.01342 

.04029 

.04029 

(.00861) 

.00017 

.00354 

(.00323) 

.00405 

.00310 
( .00350) 

22.05 

.467 

-.06672 

(-.0584) 

-.01292 

.00008 

.02162 

.02184 

(.03083) 

-.00254 

.00075 

(-.00022) 

.00458 

.00503 

(.00547) 

11.05 

.348 

-.04175 

(-.0324) 

-.00260 

-.00535 

.00138 

.00165 

(.01312) 

.00200 

.00247 

(.00158) 

.00264 

.00302 

(.00335) 

11.06 

.349 

-.05015 

(-.0324) 

.00195 

.00178 

-.01449 

-.01411 
( .00446) 

.00465 

.00438 
( .00209) 

.00241 

.00216 

(.00315) 

82.10 

.250 

.01905 

(.0146) 

- .01584 

.00124 

.06376 

.06401 

(.06005) 

.00162 

.00095 
( .00101) 

.00283 

.00308 

(.00310) 

83.08 

.248 

-.00421 

(-.0397) 

.00597 

.00501 

; .07038 

.07043 

(.03773) 

-.01320 

-.00989 

(-.00472) 

.00510 

.00553 

(.00484) 

72.06 

.472 

.03014 

(.0063) 

.01304 

.01708 

^.04121 

.04086 
( .00911) 

-.00107 

.00334 

(.00314) 

.00560 

.00346 

( .00359) 

73.08 

.468 

-.03533 

(-.0594) 

.OHIO 

-.00611 

.04991 

.04815 

(.02186) 

-.00990 

-.00372 

(.00002) 

.00809 

.00698 

(.00530) 


* Nonparenthesized values pertain to trimmed analytic resiilts; parenthesized values 
pertain to untrimmed analytic results. 




TABLE X 


SUMMARY OF FULL-SCALE DESIGN FORWARD FLIGHT TRIM CASES 


Case 

1 

2 

3 

blade configuration 

subcritical 

super c r i tic al 

subcritical 

torque-tube bending 
stiffness 

nominal. 

^ 25^ nominal 

nominal. 

radial pushrod 
attachment point, m. 

.533 

.483 

.635 

(i = 1 

.118 

-.013 

.383 

A 0 /Aq 1 = 2 

-.330 

.029 

-1.053 

(i = 3 

00 

-.024 

2.301 

^.75R> Aeg. 

12.07 

11.25 

11.82 

A^, deg. 
s 

- 1.19 

-1.89 

-1.36 

, deg. 

Q 

9.69 

8.07 

9.47 

X 

-.0609 

-.0642 

-.0581 
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TABLE XI 


AERODYNAMIC STATIC STABILITY DERIVATIVE MATRICES, [s] 


trim case 1 

(basic subcritical configuration) 



-.002090 

-.009384 

.015010 

.009280 

-.224900 

-.000060 

.000522 

.000378 

-.000432 

-.051680 

-.000177 

-.001431 

.001708 

.000685 

.000306 

.000790 

-.000089 

.000166 

.000022 

.002155 

-.000067 

-.000245 

.001318 

.000302 

-.009352 

.001148 

-.000076 

.000034 

-.000112 

.001190 __ 

trim case 2 

(basic supercritical configuration 


-.000800 

-.009808 

.016630 

.009280 

-.196000 

-.000079 

. 000444 

.000653 

-.000315 

-.046590 

-.000066 

-.001364 

.001225 

.000427 

.004696 

.000902 

-.000299 

.000857 

.000140 

-.005363 

-.000046 

-.000186 

.001176 

.000214 

-.007553 

.001137 

.000297 

-.000379 

-.000340 

.002126 _ 


C. trim case 3 (subcritical configuration with modified pitch- flat coupling) 


-.004567 

-.008454 

.014690 

.009394 

-.257600 

.000052 

.000483 

.000515 

-.000423 

-.051220 

-.000490 

-.001310 

.001583 

.000723 

-.002830 

.000594 

-.000477 

.000602 

.000236 

.000931 

-.000165 

-.000281 

.001376 

.000362 

-.015480 

.000950 

-.000538 

.000572 

.000147 

.000852 ^ 
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TABLE XII 


TRIM DERIVATIVE MATRICES, [t] 


A. trim case 1 (basic subcritical configuration) 



3-3 

-86.9 

-130.4 

1238.2 


115.^ 

747.4 

-1188.6 

96.5 


24.2 

896.5 

157.7 

167.9 


1.3 

-5.4 

-11.3 

0.8 

trim 

case 2 

(basic supercritical configuration) 



-17.5 

-449.0 

-247.7 

1036.1 


56.4 

427.0 

-1002.5 

13.4 


15.0 

787.7 

128.3 

91.5 


1.0 

-7.6 

-10.0 

-0.5 

trim 

case 3 

(subcritical configuration 

with modified pitch-flat coupling) 


6.3 

-52.6 

-532.5 

1297.5 


123.2 

570.3 

-1287.7 

-165.3 


18.2 

789.3 

113.4 

164.1 


1.3 

-5.6 

-12.0 

1.69 
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Mode 


Flexbeam Bending Stiffness 



Nominal 

(EIf)pb" 

0.85 X Nominal 

(eie)fb 

0.85 X Nomii^ 

Subcritical Design 





I.O8OP 

I.076P 

I.O8OP 

3 3 

2.17UP 

2.7OIP 

2 . 7 i 4 p 

4 . 843 P 

4 . 8 i 4 p 

4 . 834 P 


I.328P 

I.328P 

I.27OP 

tOvg • ■ . 

7.28IP 

7.28IP 

7.IIIP 

Supercritical Design 




Wwi 

I.08OP 

I.076P 

I.08OP 

a>V2 

2 . 714 p 

2.70IP 

2 . 714 p 


4 . 843 P 

4 . 8 i 4 p 

4 . 8 s 4 p 


0 . 794 p 

0 . 794 P 

O.763P 

««>V2 

4.138P 

4.138P 

4.063P 












TABLE XIV 


SELECTED EIGENSOLUTIONS FOR MODEL ROTOR WITH CANTILEVERED TORQUE TUBE 


Trailing Edge Push-Rod Attachment at 
flR = 99*1 m/sec., ^ = 0 

®.T5R “ 


= = 0, X (inflow) = 


.09R. 

-.005 


Case & Condition 

V 


Modal Content 

I. /u,= 0.1 

zero initial conditions 

-.244 

-.026 

-.210 

-.367 

-.176 

1.194 

1.356 

2.947 

3-^33 

5.794 

IF, 2T 
IE 

2T, 2F, IT 
2T, IT 
IT, 2T, 3F 

rr. /«■ = 0.1 

nominal initial conditionE 
( q = q = 0.02) 

Wj Vi 

-.227 

-.019 

+.052 

-.648 

-.181 

1.292 

1.382 

3.145 

3.130 

5.792 

IF, 2T, IE 
IE, 2T, IF 
2T, IT, IF, 2F 
2T, IT, IF, 2F 
IT, 2T, 3 P 

III. 0.47 

nominal initial conditions 

-.200 

».039 

+.054 

-.666 

-.187 

1.392 

1.393 
3.198 
3.166 
5.821 

IF, 2T, IE 
IE, IF, 2T 
2T, IT, IF, IE 
2T, IT, IF, IE 
2T, IT, IF, IE 

♦Eigenvalue, X = o' ± io) 
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TABLE XV 


EIGENSOLUnONS FOR FULL SCALE DESIGN 
CONFIGURATIONS AT TRIM CONDITIONS 


OR “ 198.1 m/sec, = 0 


(initial conditions )» (initial conditions) 




. 02 ) 


C£l80 

cr 

* 

0) 

Modal Content 

I. Subcritical 

(Configuration No.l) 

-.U99 

-.067 

-.389 

-M 3 

-.1U7 

-2.713 

1.201 

1.403 

2.836 

4.789 

7.575 

10.059 

IF, IE, IT 
IE, IF, IT 
2F, IT, IE, IP 
3F, IE, IT, 2P 
IT, 2E, IP 
IT 

II. Supercritical 

(Configuration No. 2) 

- .058 
-.563 

-.kl 3 

-.081 

-.505 

-2.646 

.847 
1.304 
2.808 
4 .078 
4.944 
10.081 

IE, IF 

IF, IE 

2F, IE, IP, IT 
IT, 2E, IP, 3F 
3P, IE, IT, 2E 
IT 

♦Eigenvalue X=o’±1cd 
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YOUNG’S MODULUS. E . GN/m 


^ 60 % FILL 



j Figure 1. -Typical Combinations of E and G. 

j 


! 


I 





'S 



Figure 2. - Composite Bearingless Rotor 


Figure 3 • ” Fatigue Test Setup 










ALTERNATING SHEAR STRESS BASED ON 
INITIAL MODULUS, MN/m^ 


Figure 5 . - Torsion Fatigue Results . 
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+ 52.7 MN/m^ PURE TORSION 
5 X lo‘ CYCLES 




Figure 6. - Failed Graphite/Epoxy Fatigue Specimen 


















BLADE THICKNESS, cm 



BLADE CHORD, cm 


Figure 10. -Typical Cross Section of Model Blade 







Pinned- Pinned Torque Tube. 
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0.5 


- 1.0 


Figure lU. -Calcu! 


EDGEWISE & TORSION MODE SHAPES, & 70 FLATWISE MODE SHAPES, 7 




BLADE NATURAL FREQUENCIES, co, Hz 



TIP SPEED, m/sec 


Figure 15. - Variation of Blade Natural Frequencies with Tip Speed, 
Pinned- Pinned Torque Tube. 
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300i 





THRUST COEFFICIENT/SOLIDITY, Cj/or 





rORQUE COEFFICIENT/SOLIDITY, ■ Cq/o 


Figiire l8. -Model Rotor Hover Perfonnance, Cantilevered Torque Tube, 
Trailing Edge Push- Rod Attachment at 0.06 R, cr = 0.04o8. 







THRUST COEFFICIENT/SOLIDITY, Cj/O 








rr COEFFICIEI 







THRUST COEFFICIENT/SOLIDITY, Cj/a 







THRUST COEFFICIENT/SOLIDITY, Cjh 
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S 



Figure 22. -Effect of Blade Pitch Setting and Push Rod Attachment Point 
Location on Model Rotor Hovering Thrust, flR = 99 m/sec, a- = 0.o4o8. 




TORQUE COEFFICIENT/SOLIDITY, Cq/o DRAG COEFFICIENT/SOLIDITY, Cq/jt lIFT COEFFICIENT/SOLIDITY, Cja 



Figure 23. -Effects of Input Collective Angle and Shaft Angle Variations on 
Forward Plight Performance of MDdel Rotor with Pinned-Pinned 
Torque Tube, fi = 0.25, cr = O.oSl?. 
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TORQUE COEFFICIENT/SOLIDITY. Cq/o DRAG COEFFICIENT/SOLIDITY, Cp/a LIFT COEFFICIENT/SOLIDITY, CJa 





TORQUE COEFFICIENTySOLIOITY. Cq/c COEFFICIENT/SOLIOITY. Cp/a LIFT COEFFICIENT/SOLIOITY. Cl/« 





SHAFT ANGLE, os, deg 


Figure 2$. -Effects of Input Collective Angle and Shaft Angle Variations on 
Forward Fli^t Performance of Model Rotor with Pinned- Pinned 
Torque Tube, fjL = 0.^7, <r = 0.08l7- 
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TORQUE COEFFICIENT/SOLIDITY, Cq/o DRAG COEFFICIENT/SOLIDITY, Cq/o LIFT COEFFICIENT/SOLIDITY, Cl/o 





SHAFT ANGLE, 03, deg 

Figure 2 ?. -Effects of Input Collective Angle and Shaft Angle Variations on 
Forvreird Flight Performance of Model Rotor with Csmtilevered ; 
Torque Tube, Leading Edge Push- Rod Attachment @ O.O7 R, 
jx = 0.25, cr = O.0U08. 
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1^2 


SHAFT ANGLE, as, deg 

Figure 28. -Effects of Input Collective Angle and Shaft Angle Variations on 
Forward Plight Performance of Model Rotor with Cantilevered 
Torque Tube,iTrailing Edge Push- Rod Attachment @ 0.08 R, 

/A * 0.25, a- - 0.04o8. 




TORQUE COEFFICIENT/SOLIDITY, Cq/0 DRAG COEFFICIENT/SOLIDITY.Cq/w LiFTCOEFFICIENT/SOUDITY.CL/a 



-12 -10 -8 - 6 - 4-2 0 +2 

SHAFT ANGLE, as, deg 

Figure 29. -Effects of IhTait Collective y^le and Shaft )higle Variations on 
Forward Flight Performance of Model Rotor with Cantilevered 
Torque Tube » Leading Edge Push- Rod Attachment @ 0.09 
IX - 0.25, cr = 0.0408. 
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TORQUE COEFFICIENT/SOLIDITY, CQ/a DRAG COEFFICIENT/SOLIDITY.CQ/a LIFT COE FF I Cl ENT/SOLI DITY,C|_/a 


5 ] 


0 = 11 . 2 ° 


SHAFT ANGLE, as, deg 

^ Effects of Input Collective Angle and Shaft Angle Vaxlatlons on 
Forward Flight Performance of Model Rotor with Cantilevered. 
Torque Tube, Leading Edge Push- Rod Attachment @ O.09 R, 

M = 0.35, o- = 0.0408. 






TORQUE COEFFICIENT/SOLIDITY, Cq/w DRAG COEFFICIENT/SOLIDITY. Cp/ a upT COEFFICIENT/SOLIDITY, C^l a 





Figure 32. -Effects of Input Collective Angle and Shaft Angle Variations on 
Forward Flight Performance of Model Rotor with- Cantilevered 
Torque Tube, Trailing Edge Push- Rod Attachment @ 0.06 R, 

IX = 0.47, o- = 0.0408. 
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TORQUE COEFFICIENT/SOLIDITY, Cq/o DRAG COEFFICIENT/SOLIDITY^CQ/a LIFT COEFFICIENT/SOLIDITY, ^\^fa 





SHAFT ANGLE, as, deg 

Figure 33. -Effects of Input Collective Angle and Shaft Angle Variations on 
Forward Flight Performance of Model Rotor with Cantilevered 
Torque Tube » Leading Edge Push- Rod Attachment @ O.09 R» 

IX = 0.47, a- = O.0U08. 













1/2 PTP FLATWISE BENDING MOMENT, Mp. Nm 



Figure 36. - Effects of Advance Ratio and Blade Loading on Vibratory Outboard 
Flatwise Bending Moment of Model Rotor with Pinned- Pinned Torque 
Tube , r = .42 R, ^ - O.O817. 
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0 

o 

0.25 

13.7° 

a 

0.35 

13.6° 

A 

0.47 

13.6° 


0.25 

7.5° 

o' 

0.35 

7.5° 

A* 

0.47 

7.3° 



Figure 37 • -Effects of Advance Ratio and Blade Loading on Vibratory Aboard 
Flatvdse Stresses of Model Rotor with Cantilevered Torque Tube, 
Trailing Edge Push-Rod Attachment @ 0.06 R, r = 0.06 R, <r = 0.0408. 
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Figure 38. -Effects of Advance Ratio and Blade Loading on Vibratory Inboard 
Edgewise Stresses of Model Rotor with Cantilevered Torque Tube, 
Trailing Edge Push-Rod Attachment @ 0.06 R, r = 0.06 R, cr= 0.0^ 





1/2 PTP FLATWISE BENDING MOMENT. Mp , Nm 



THRUST COEFFICIENT/SOLIDITY, Cj/a 


Figure 39. -Effects of Advance Ratio and Blade Loading on Vibratory 
Outboard Flatwise Bending Moment of Model Rotor with 
Cantilevered Torque Tube Trailing Edge Push-Rod Attachment 
@ 0.06 R, r = 0.42 R, <r= 0.o4o8. 
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9 

o 

0.28 

11.2® !. 

D 

0.36 

11.2® ■ 

A 

0.47 

\ 11.2®.'V.' 

cr 

0.25 
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0.35 

13.4® ■ 

a 

0.47 

8.S® « ■ 
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0.47 
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t V-/ 


6 — ^ ^ ySi 
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2 — POINTS WITH ARROWS DENOTE VALUES FOR 
I ALTERNATE REVOLUTIONS |Of CONDITIONS 

I EXPERIENCING 1/2 P OSCILLATIONS 


0.02 0.04 0.08 

THRUST COEFFiCiENT/SOLIOiTY, Cj/o 


•Effects of Advance Ratio and Blaite Loading on Vibratory lUboard 
F^twlse Stresses of Model Rotor with Cantilevarad Torque Tuiba» 
leading Edge Push-Bod Attachnant 9 0.09 R* r ■ 0.06 H, <r « O.OboS 



^ 1/2 PTTEfXIi-Wlsi STRESS, MN/m^ 


14 


12 
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THRUST C(XFFICIENT/SOLIOITY. Cj/a 

, .PiguTtt M* ’-EffRct* of Advrac* Ratio and Blada Loading, on Vibratozy 

' Babbard EdgaidM Straaaas of Nodal Rotor with Caatilavarad 
■ ' " *f6r^'!Pid)a, Ti^ Bii^' PuahiRbd AttachB«nt;it 0.06 R, 

■ ./■ - r « 0.06 R, am O.04o8.' ■ ' -H. 
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1/2 FTP FLATWISE BENDING MOMENT.Mf;, Nm 



Figure k2. -Effects of Advance Ratio and Blade Loading on Vibratory 
Outboard Flatwise Bending Mament of Model Rotor with 
Cantilevered Torque Tube Trailing Edge Push-Rod 
Attachment @ 0.06 R, r = 0.42 R,a-=0.0408. 
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EDGEWISE STRESS, Og, MN/m^ FLATWISE STRESS, Op, MN/m^ 




(a) SHAFT ANGLE PERTURBATIONS 


Figure 43. -Perbiirbations of Azimuthal Variations of Inboard Flatwiae 
and Edgewise Stresses for Model Rotor with Finned-Finned 
Torque Tube, /x = 0.25, or= 0.0408. 
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AZIMUtH POSITION, tfi. dq 
(b) COLLECTIVE PITCH PERTURBATIONS 

Figure U3> ** (Continued) - 





EDGEWISE STRESS^ MN/m^ fUtwISE S1RESS. Of. MN/m^ 












azimuth position, )/<, deg 



(a) SHAFT ANGLE PERTURBATIONS 

Figure UU. -Perturbations of Azimuthal Variations of Inboard Flatwise and 
Edgewise Stresses for Model Rotor with Pinned-Pinned 
Torque Tube, ix = 0.35» cr = 0.0408. 
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AZIMUTH POSITION, deg 


(b) COLLECTIVE PITCH PERTURBATIONS 


Figtire hh. - (Continued) 






EDGEWISE STRESS, a^, MN/m^ FLATWISE STRESS, Op, MnW 




(c) LONGITUDINAL CYCLIC PERTURBATIONS 


Figure hk. - (Continued) 




AZIMUTH POSITION, deg 


(d) LATERAL CYCLIC PERTURBATIONS 


Fig\ire — (Concluded) 





EDGEWISE STRESS, ag, MN/m2 FLATWISE STRESS, op. MN/m^ 




AZIMUTH POSITION, deg 
(a) SHAFT ANGLE PERTURBATIONS 


Flgtire 45* -Perturbations of Azimuthal Variations of Inboard Flatwise 
and Edgewise Stresses for Model Rotor with Cantilevered 
Torque Tube, Trailing Edge Push- Rod Attachment @ 0.06 
R, /u = 0.25, o- = 0.0408. 
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AZIMUTH TOSITION, deg 
(b) COLLECTIVE PITCH PERTURBATIONS 


Figure U5. - (Continued) 
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EDGEWISE STRESS, a^. MN/m^ FLATWISE STRESS, ap, MN/m^ 












(a) SHAFT ANGLE PERTURBATIONS 

Figtire h6. -Perturbations of Azimuthal Variations of Inboard Flatwise and 

Edgewise Stresses for Model Rotor with Cantilevered Torque Tube, 
Trailing Edge Push-Rod Attachment @ 0.06 R, /4 = 0.35» = 0.04o8. 
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EDGEWISE STRESS, aE,MN/m2 FLATWISE STRESS, Op, MN/m^ 














0 90 180 270 360 

AZIMUTH POSITION,!;^, deg 


(a) SHAFT ANGLE PERTURBATIONS 


Pig\ire U7. -Perturbations of Azimuthal Variations of Inboard Flatwise and 

Edgewise Stresses for Madel Rotor with Cantilevered Torque- Tube, 
Leading Edge Push- Rod Attachment at 0.09 E, = 0.3^, cr = 0.0408. 
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AZIMUTH POSITION, deg 


(b) COLLECTIVE PITCH PERTURBATIONS 


Figure UT. - (Continued) 




















LIFT COEFFICIENT, Cj 










EDGEWISE BENDING STRESS, OE,MN/m2 



Figure 50 . -Experimental aind Analytic Time-Histories of Inboenrd Edgewise Bending 
Stresses, r = 0.06 R, Case 25 * 03 , C^/<r = 0 . 0783 > /x = 0 . 25 , b = 2 . 






FLATWISE BENDING STRESS, a p, MN/m^ 



Figure 52. -Experimental and Analytic Time- Histories of Inboard Flatwise Bending 
Stresses, r = 0.06 R, Case 87 . 08 , C /cr = 0.07^1, fi= 0.25, b = 4. 
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FLATWISE BENDING MOMENT/ Nm 


0.02 



Esgperimezrtal and Analytic Time-HLstories of Outboard Flatwise Bending 
Ncnents, r s 0.42 R, Case 87 . 08, C./o* = 0.0741, a s 0.2^, b s 4. 


Figure 54 . - 
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BLADE AZIMUTH POSITION, 4i, DEG 


Figure 55. -Experimental and Analytic Time-Histories of Inboard Flatwise Bending 

Stresses, r = 0.06 R, Case 22.03, C /cr = O.0708, u, = 0.25. 
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EDGEWISE BENDING STRESS,, ag, MN/m2 



Figure 56 . -Experimental and Analytic Time-Histories of Inboard Edgewise Bending 
Stresses, r = O.O 6 R, Case 22.03, C /<r = 0,0708, fi = 0.25. 
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FLATWISE BENDI 



BLADE AZIMUTH POSITION, if/, DEG 


Figure 57. -Experimental and Analytic Time-Histories of Outboard Flatwise Moments, 
r = 0.42 R, Case 22.03, OjJa = O.0708, p, = 0.25. 






FLATWISE BENDING STRESS, o p.MN/m^ 


2 



Figure 58. -Experimental and Analytic Time-Histories of Inboard Flatvise Bending 
Stresses, r = O.06 R, Case 22.04, CJ<r = O.0503, /*= 0.35« 
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0 90 180 270 360 


BLADE AZIMUTH POSITION. DEG 


Figure 62. - Essperimental and Analytic Time-HLstories of Inboard Edgewise Bending 
Stresses, r = 0.06 R, Case 21.05, C^/<r = 0.0403, fi= 0,47. 



FLATWISE BENDING MOMENT, Nm 



Figure 63. -Experimental and Analytic Time-Histories of Outboard Flatwise Bending 

Moments, r = 0.42 R, Case 21.05, C /<r = 0.0403, yt, = 0.47. 
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Figure 6k. -Experimental and Analytic Time-Histories of Inboard Flatwise Bending 
Stresses, r = 0.06 R, Case 22.05, = 0.0216, fx = 0.4y. 


H 


EDGEWISE BENDING STRESS, a^.lMN/m^ 
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BLADE AZIMUTH POSITION. \l/. DEG 


Figure 65 . -Experimental and Analytic Time-Histories of Inboard Edgewise Bending 
Stresses, r = O.O 6 R, Case 22.05, = 0.0216, /j, = 0.47. 



FLATWISE BENDING MOMENT, Nm 


0.02 



Figure 66. -E:^rimental and Aneilytic Time-Mstories of Outboco^ Flatvdse Bending 
t Manents, r = 0.42 R, Case 22.05, C_/<r = 0.0216, u= 0.4?. 
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BLADE AZIMUTH POSITION, DEG 


Figiire 68. -Experimental and Analytic Time-Histories of Inboard Edgewise Bending 

Stresses, r = 0.06 R, Case 11.05, C_/cr = 0.0014, fj, = 0.35. 
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BLADE AZIMUTH POSITION, DEG 


Figure 70. -Esgperlmental and i\nalytic Time- Histories of Inbosird 
Edgewise Bending Stresses, r = 0.06 R. Case 11. 06, 






Figure 71* -Experimental and Analytic Time-Histories of Inboard Flatwise Bending 
Stresses, r = 0.06 R, Case 82.10, C /cr = O.0638, a = O.25. 
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BLADE AZIMUTH POSITION, DEG 


Figure 72. -Experimental and Analytic Time-Histories of Inboard Edgewise Bending 
Stressei^, r = 0.06 R, Case 82.10, = O.0638, jj,= O.25. 









BLADE AZIMUTH POSITION, ii. DEG 


Figure TU. -Experimental and Analytic Tiiue-Histories of Flexbeam Torsional Deflections, 

Case 82.10, C,/<r ■ O.O 638 , fi & 0.2$. 
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FLATWISE BENDING STRESS, a p, MN/m2 



BLADE AZIMUTH POSITION, DEG 

Figure 75 • - £3g;>exl]nental and Analytic Time-IEstories of Iaiboa3rd Flatwise Bending 
Stz«sses, r = 0.06 R, Case 83 . 08, Cj/<r * 0.0704, ft = 0.25. 



179 



Figure 76. -Experimental and Analytic Time-Histories of Inboard Edgewise Bending 
Stresses, r = O.06 R, Case 83.08, = 0.0704, fi = 0.25. 








TORSIONAL DEFLECTION. DEG 



BLADE AZIMUTH POSITION, DEG 

Figure 78. -Ejgperinental and Analytic Ume-HLstories of Flexbeam Torsional Deflections, 
Case 83.08, Cj^/<r = 0 . 0704 , ^ = 0.25. 
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BLADE AZIMUTH POSITION, DEG 


Figure 79* - 


Experimental and Analytic Time-Histories of Inboard Flatvd.se Bending 
Stresses, r = 0.06 R, Case 72-06, C^/cr = 0.04l2, fi = 0.47. 






Figure 80. - Eacperimental and Analytic Time-Histories of Inboard Edgewise Bending 
Stresses, r = 0.06 R, Case 72.06, C /c = 0.04l2, /x= 0.47« 
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l-LATWISE BENDING MOMENT, Nm 



BLADE AZIMUTH POSITION, DEG 

Figure 8l. - Experimental and Analytic Time-Histories of Outboard Flatwise Bending 
Monents, r = 0.U2 R, Case 72.06, C_/cr = 0.04l2, /u. = 0.47. 




TORSIONAL DEFLECTION, DEG 



Figure Q2. - Esgperimental and ilnalytic Time- Histories of Flexbeam TorslonaJ. 
Deflections, Case 72.06, C^/V = 0,04l2, fi = 0.47. 
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Figiire 03. -ExperimentaJ. and Analytic Time-Ifi.stories of Inboard Flatvdse Bending 
Stresses, r = 0.06 R, Case 73.08, C /<r = 0.04l2, /x = 0.4?. 
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Figure 8i^. -Eatperimental and Analytic Time-Histories of Inboard Edgewise Bending 

Stresses, r = 0.06 R, Case 73*08, 0 / or = 0.04l2, /x = 0.47. 
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BLADE MASS AND INERTIA DISTRIBUTIONS 



Figure 88. -Full Scale Weight and Inertia Distributions 
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Flg^lre 90. -Full Scale Utorsion Stiffness Distributions. 
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TORQUE TUBE SECTION PARAMETER DISTRIBUTIONS 
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Figure 91* -Full Scale Torque Tu*be Characteristics, Cantilevered and 
Snuhber Designs 
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Figure 95* “Full Scale Blade Uncoupled Frequencies of Low NumUered Modes 
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Figure 97* —Variation of Coupled Frequencies with Tip Speed, 
Subcrltlcal Design 
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Figure 100. -Variation of Coiipled Frequencies with Tip Sx>eed 
Supercritical Design 




BLADE COUPLED NATURAL FREQUENCY, «, Hz 



Figure 98. - Variation of Coupled Prequenclea with Collective Pitch 
Angle, Suhcrltlcal Design 
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Figure 102. -Variation of Coupled Frequencies with Torque Tube 
Bending Stiffnesses^ Supercritical Design 
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Flgtire 103 . — Spanwlse Variation of Coupled Mode Shape Components 

0 , Subcritical Design 
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NONVACUUM COUPLED MODE CRITICAL DAMPING RATIO, 
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COLLECTIVE PITCH ANGLE, 0 75 , deg 

Figure 107. —Variation of Nonvacuum Coupled Mode Damping with 
Collective Angle, Nonstalled Aerodynamics, 
Subcritical Design 





NONVACUUM COUPLED MODE CRITICAL DAMPING R A 










.8 


li= 0 (HOVER) 


nR = 198 m/$ec 
= 1.33 

NOMINAL TORQUE TUBE 
BENDING STIFFNESS 


- STATIC STALL 
ANGLE OF 11.3° 
ENCOUNTERED 
AT 75% SPAN 




STABLE 


u 


UNSTABLE 


(IF, IE) 


(2F. IE) 


(IT, IE. 2E) 


(IE. 3F) 


(IF, IE) 


COLLECTIVE PITCH ANGLE, 6 75, deg 

Figure 109 • - Varliitlon of Noirvacuum Couj>led Mode Damping with 
Collective Angle^ Including Stall Characteristics, Suhcrltlcal. Design 
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Pigixre 110. - Variation of Nonvacuum Coupled Mode Damping with 

Collective Jtagle, Including Stall Characteristics, Supercritical Design 

















NONDIMENSIONAL 


Figure lilt, - Spanwise Variation of 



Flatwise Stresses, Trim Cases 1 and 2 








EDGEWISE BENDING STRESS, o^, MN/m2 



Fig\ore 115 • ~ ^anwise Variation of FlexTjeam Edgewise 
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Figvire ll6. - Azimuthal Variation of Critical Bending Stresses, Trim Cases 1 and 2 
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Figure 119a* - Effects of Complete Modelling of Structural 
Redundancy on Stresses 
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Figure 120. - Incremental Ramp- Step and Pulse Input Control Angle Time-Histories 
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FlgiJre 121a 


Transient Bending Stresses due to Collective Angle 
Inputs, Subcrltlcal Design 
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Figiire 1221). “ Transient Flexbeam Torsion Stresses and Push-Rod 

Loads due to Collective Angle Inputs, Supercritical 
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Figure 12313. - Transient Flexbeam Torsion Stresses and Push-Rod 
Loads due to Longitudinal Cyclic Angle Inputs, 
Suhcritical Design 
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Flgvire 12lrt). - Transient Flexbeam Torsion Stresses and Push-Rod 
Loads due to Longitudinal Cyclic Angle Inputs, 
Supercritical Design 
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Figure 12613. - Transient Flextieam Torsion Stresses and Push-Rod Loads 

due to Lateral Cyclic Angle Inputs, Supercritical Design 
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Figure 130. - Effects of Unsteady Airloads Formulation on Variations 

of Vilsratoiy Bending Stresses and Push-Rod- Load with Shaft Angle 
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MODAL EQUIVALENT CHARACTERISTIC NONDIMENSIONAL IDENTIFIED MODAL 
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Figure 131. - Variation of Identified Modal Frequencies and Characteristic 
Exponents with Advance Ratio, Cantilevered Torque Tuhe 
Model Rotor, Trailing Edge Push-Rod Attachment at .080R 
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Figure 133. - Variation of Frozen Azimuth Eigenvalues with Tip Speed, 
Cantilevered Torque Tuhe Model Rotor, Trailing Edge 
Push-Rod Attachment at . 080 R 
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Figure 13*** - Variation of Identified Modal Frequencies and Characteristic 
Exponents with Tip JSpeed, Cantilevered Torque Tube 
Model Rotor, Trailing Edge Push-Rod Attachment at .080R 
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OPEN SYMBOLS: SUBCRITICAL CONFIGURATION 
CLOSED SYMBOLS: SUPERCRIT: CONFIGURATION 


F: FLAPPING RESPONSE AT TIP 
I; INPLANE RESPONSE AT TIP 
T: TORSION RESPONSE AT TIP 
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Figure 136. -Variations of Identified Modal Frequencies and Characteristic Exponents 

with Mass Ratio and Shaft Angle > Suhcrltical and Supercritical Full Scale Designs 
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OPEN SYMBOLS: SUBCRITICAL CONFIGURATION 
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Figure 138. - Variation of Identified Modal Frequencies and Characteristic 
Exponents with Push-Rod Attachment Location, Suhcritical and 
Supercritical Full Scale Designs 
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